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ON THE RATE OF CONVERGENCE TO STATION ARITY OF THE M/M/N 
QUEUE IN THE HALFIN-WHITT REGIME 

By David Gamarnik and David Goldberg 
MIT 

We study the rate of convergence to stationarity of the M/M/n 
queue in the Halfin-Whitt regime. We prove that there is an inter- 
esting phase transition in the system's behavior, occurring when the 
, excess parameter B reaches B* « 1.85772. For B < B* , the expo- 

nential rate of convergence is above B* it is the solution to an 
}—{ , equation involving the parabolic cylinder functions. We also bound 

the prefactor governing the rate of convergence uniformly over n when 
B < B", and use our bounds to derive a rule-of-thumb for determin- 
ing the time it takes a severely overloaded (underloaded) queueing 
system to return (probabilistically) to the steady-state. 

I 1. Introduction. Parallel server queueing systems can operate in a variety of regimes that 

balance between efficiency and quality of offered service. This is captured by the Halfin-Whitt 
• (H-W) regime, which can be described as critical with respect to the probability that an arriving 

, customer has to wait for service. Namely, in this regime the stationary probability of wait is bounded 

away from both and 1 as the number of servers grows. Although studied originally by Erlang [16] 
and Jagerman [25], the regime was formally introduced in [22] by Halfin and Whitt, who established 
the criticality of the probability that a customer waits in terms of the square-root spare capacity 
^ I rule for large multiserver queues with exponentially distributed service times. They also showed 

that as the size of the multiserver queue grows, the queueing process (properly scaled) converges 
. weakly to a non-trivial positive recurrent diffusion (the H-W diffusion). Similar weak-convergence 

I results were subsequently obtained for more general multiserver systems [34], [26], [32], [19], with 

' the most general results appearing in [36]. This line of research has been motivated by questions 

related to the management of call centers, since it is generally accepted that well-run call centers 
operate in the H-W regime [1],[20]. 

Since the steady-state behavior of the M/M/n queue in the H-W regime is quite simple in 
practice [22], while the transient dynamics are more complicated [22], it is common to use the 
steady-state approximation to the transient distribution. Thus it is important to understand the 
\ quality of the steady-state approximation. The only work along these lines seems to be the recent 

paper [30], which studies the transform of the H-W diffusion and proves several results analogous 
to our own for the H-W diffusion. The key difference is that in this paper we study the pre-limit 
M/M/n queue, not the limiting diffusion. We note that the relevant transform functions were also 
studied in [2], although in a different context. 

The question of how quickly the positive recurrent M/M/n queue approaches stationarity has a 
rich history in the queueing literature. In [33], Moorse derives an explicit solution for the transient 
M/M/1 queue, and discusses implications for the exponential rate of convergence to stationarity. 
Similar analyses are carried out in [10] and [37]. Around the same time, two different research 
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groups [31], [27] worked out powerful and elegant theories that could be used to give the tran- 
sient distributions for large classes of birth-death processes (b-d-p). The transient probabilities 
are expressed as integrals against a spectral measure (pix) that is intimately related to the eigen- 
values of the generator of the b-d-p. Karlin and McGregor (K-M) devote an entire paper [28] 
to the application of their theory to the M/M/n queue. For a positive recurrent b-d-p Z{t)^ let 



Pij{t) = Pr(Z(t) = j|Z(0) = i), P{t) = {Pij;i,j G Z+), and Pj(oo) = liuit^^ Pij{t) denote the 



steady-state distribution. Then K-M proved 

Theorem 1. [28] For any positive recurrent M/M/n or M/M /oo queue, there exist a set of 
polynomials {Qi{x)}, a set of strictly positive constants {cj}, and a measure (j){x) with non-negative 
support such that 



K-M also comment on the relationship between the rate of convergence to stationarity and 
the support of (pix) [28]. This relationship was later made rigorous and expanded on in a series 
of papers [4], [5], [40], [6]. Recall that the spectral gap 7 of a b-d-p is the absolute value of the 
supremum of the set of strictly negative real eigenvalues of the associated generator matrix A, 
where ^P{t) = AP{t) [17], [6]. If no such eigenvalues exist, we set 7 = 00. For a complex-valued 
function /(x), we let Z (f) (^Z~^ (f)^ denote the infimum of the set of (strictly positive) real zeros of 
f{x). Set Z(/) = 00 if f{x) has no (strictly positive) real zeros. We then have 

Theorem 2. [28], [4], [5], [40], [6] For any positive recurrent M/M/n or M/M/00 queue, 7 e 
(0, 00). For all i and j, limj^oo —t~^ log \Pij{t) — -Pj(oo)| > 7. For at least one pair of 
limt_!.oo —t~^ log \Pij{t) — Pj(oo)| = 7. Furthermore, 7 = Z~^{cl)). 

It is well-known that for the positive recurrent M/M /I and M/M/ 00 queues, 7 can be computed 
explicitly. In particular. 

Theorem 3. [28] For the positive recurrent M/M/1 queue with arrival rate A and service rate 
7 = (\/A — ^/Jxf' , and the spectral measure 4){x) consists of a jump at 0, and an absolutely 
continuous measure on [(\/A — y^)^, (\/A + For the M/M /oo queue with arrival rate A and 

service rate 7 = /i, and the spectral measure 0(x) consists of a countably infinite number of 
jumps, at the points {kfi;k € Z^}. 

Unfortunately, for the general positive recurrent M/M/n queue, the known characterizations 
for 7 are cumbersome and hard to use. The difficulty is that for the positive recurrent M/M/n 
queue with arrival rate A and service rate the spectral measure 4>(x) consists of three parts [28]. 
The first part is a jump at 0, which corresponds to the steady-state distribution [28]. The second 
component is an absolutely continuous measure on the interval [(\/A — ^n/x)^, (-v/A + y^^)^], whose 
density is described in [28]. The third component consists of a set of at most n (but possibly zero) 
jumps, which all exist on (O, (\/A — -^/n/I)^) [28]. The complexity of determining 7 arises from the 
difficulty of locating these jumps [40]. In [28], the set of jumps is expressed in terms of the zeros of 
a complicated polynomial equation. 

Significant progress towards understanding these jumps was made in a series of papers by van 
Doorn [38], [39], [40], [41]. Van Doorn used the K-M representation and the theory of orthogonal 
polynomials to give several alternate characterizations and bounds for the spectral gap of a b-d-p. 
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and applied these to the M/M/n queue. He also showed that for each fixed n there is a transition 
in the nature of the spectral measure of the M/M/n queue as one varies the traffic intensity [38]. 



Theorem 4. [38] For all n > I, there exists /)* G [0,1) such that for any M/M/n queue 
satisfying > 7 = (\/A - ^/nfif' ; and for any M/M/n queue satisfying ^ < p*, 7 < 



Unfortunately, all of the characterizations (including that of p* ) given by van Doorn are again 
fairly complicated, and van Doorn himself comments in [40] that one is generally better off using 
the approximations that he gives in the same paper. Van Doom's work was later extended in [29], 
and similar results were achieved using different techniques in [46]. It was also shown in [46] that 



There are also some results in the literature for explicitly bounding the distance to stationarity 
(as opposed to just identifying the exponential rate of convergence) [46], [7]. These bounds are 
generally given in terms of an explicit prefactor attached to an exponentially decaying term (often 
containing the spectral gap). However, these bounds are generally not studied in the context of 
weak convergence, and thus may not scale desirably with n. In Section 8 we will show that two 
such bounds appearing in the literature ([46] and [7]) do not scale properly in the H-W regime. We 
note that the complexity of bounding the distance to stationarity uniformly for a sequence of b-d-p 
is related to the cutoff phenomenon for Markov chains [11], which has been studied in the context 
of queueing systems [18]. We also note that the asymptotic properties of the bounds appearing 
in [46] have recently been analyzed for several models different from those studied in this paper 



1.1. Outline of paper. In Section 2, we introduce some notation and state our main results. In 
Section 3, we prove some technical preliminaries, cite several results from [29] and [38], and prove 
a new characterization for 7^ (Lemma 3). Sections 4-7 are devoted to studying the asymptotic 
properties of Lemma 3, with each section analyzing a different term appearing in Lemma 3. This 
culminates with Section 7, in which the asymptotic results of Sections 4-6 are combined with the 
characterization for 7„ given in Lemma 3 to prove our main results (Theorem 5, Corollary 1, and 
Theorem 6). In Section 8 we prove our explicit bounds on the distance to stationarity (Theorem 
7), and compare to other bounds from the literature. In Section 9 we summarize the main results 
of the paper and present ideas for future research. We note that several of our proofs necessitate 
a considerable amount of simple algebra, and we have generally chosen to include the details, with 
an eye towards making the proofs easier to understand in their totality. 

2. Main Results. 

2.1. Definitions and notations. We now define several important quantities for the M/M/n 

queue with arrival rate A„ = n — B^/n and service rate = 1. We assume that > 0, and n is large 
enough to ensure that > and n > A„ + 1. Let Q{t) denote the queue length at time namely 

the total number of customers including those in service. Let -Pj" (i) = Pr(Q(t) = j|Q(0) = i), 

P,-(oo) ^ hm^^oo Pr,,{t). Pr,<,it) = Ei=o and P^^.(oo) ^ Ei=o Pkioo)- Let 7„ denote the 

spectral gap of the associated Markov chain. We let Hi. = X]j=i j denote the kth. harmonic number, 

and = lim„_!.oo {Hn — log(n)) ~ .5772 denote the Euler-Mascheroni constant. All logarithms will 
be base e. Unless otherwise stated, all functions are defined only for real values of x. All empty 




[42],[43]. 
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products are assumed to be equal to unity, and all empty summations are assumed to be equal to 
zero. 



2.2. The parabolic cylinder functions. We now briefly review the two-parameter function com- 
monly referred to as the parabolic cylinder function Dx{z), since we will need these functions for 
the statement (and proofs) of our results. For excellent references on these functions, see [21] Sec- 
tion 8.31 and Section 9.24, [3] Sections 3.3-3.5, and [15] Chapter 8. Let T{x) denote the Gamma 
function (see [23], Chapter 8.8). It is stated in [3] that for real x and z, Dx{z) G M, and 



(2) D,{z) 



2 2 

/o°° e-VCos(f X - zt)t''dt if x > 0; 



^ /o" e-V--*r if X < 0. 
Dx{z) takes on a simpler form for x £ Z. In particular, it is stated in [21] that for z G M, 

(3) D^i{z) = V2e^ e'^ dt, Do{z) = e~^, and Di{z) = ze~^ . 

Note that since r(— x) G (0, oo) for x < 0, (2) and (3) imply that 

(4) for all z£R and x < 0, Dx{z) > 0. 

The parabolic cylinder functions arise in several contexts associated with the limits of queueing 
models, such as the Ornstein-Uhlenbeck limit of the appropriately scaled infinite server queue [24] 
and various limits associated with the Erlang loss model [44]. We note that the parabolic cylinder 
functions have been studied as the limits of certain polynomials under the H-W scaling, using tools 
from the theory of differential equations [12], [13], [14], [35]. 

2.3. Main results. We now state our main results. We begin by describing the asymptotic phase 
transition that occurs w.r.t. the spectral gap 7„ of the M/M/n queue in the H-W regime. Let 



(5) CiB) = 



oo otherwise; 



and 

(6) ^oo{x) 

I oo 

Note that C(-B) = ^'oo(4^)- We include a plot of C(5) 



A I Dt[~(-l) +hiB + {B^- 4x)^) if D,^,{-B) / 0, X < f ; 
oo otherwise. 
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Fig 1. C(B) 



Let B* = Z+iC). Then 

Proposition 1. B* e [V2,2), and Z+(^oo) < mm(l, ^) for B > B* . Numerically, B* w 
1.85772. 

Our main result is that 

Theorem 5. The limit hm„_>.oo 7n = 7b exists for all B > 0. For < B < B* , 75 = For 
B > B*, jB = ^+(^00). 



We include a plot of 75. 



7b 




B 

Fig 2. 7B 



We note that due to the highly non-linear manner in which the steady-state probability of wait scales 
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in the H-W regime (see [22]), the case < B < B* actually encompasses most scenarios of interest. 
Indeed, it is proven in [22] that lim„_^oo -f<:„(oo) = 1 — (l + B exp{^B'^) J^^exp{—^z'^)dz)~ . As 
this limit is monotone in B, the case < B < B* includes all scenarios for which the steady-state 

probability of wait lies in the set (^{l + B* exp{^{B*)^) /^^ exp(-iz2)(iz)"\ 1^ (.038, 1). 

The following corollary may be interpreted as an asymptotic version of Theorem 4, unifying several 
earlier characterizations for the spectral gap of the M/M/n queue [28], [38], [40], [46], [29]. 

Corollary 1. The p* parameter of Theorem 4 satisfies 

lim V^(l -p^) = B*. 

n— ^oo 

The following theorem lists some properties of jb- 

Theorem 6. jb is a continuous, strictly ifiCTCdsifig fuTictiofi of B on (0, oo). lini^^QQ 
and lims^oo ^"^^jg^'^'^^ = ~h 

We now give an interpretation of Theorem 5 and Corollary 1. The M/M/n queue behaves like an 
M/M/1 queue when all severs are busy, and an M/M /oo queue when at least one server is idle. 
The phase transition of Theorem 5 formalizes this relationship in a new way. For < B < B* , 
the K-M spectral measure of the M/M/n queue in the H-W regime has no jumps away from the 
origin, and has spectral gap equal to (\/A^ — \/n)^, two properties shared by the associated M/M/1 
queue (see Theorem 3). For B > B*, the K-M spectral measure has at least one jump away from 
the origin, like the associated M/M/oo queue (whose spectral measure has only jumps and spectral 
gap equal to 1, see Theorem 3). Another interpretation is that the M/M/n queue cannot approach 
stationarity faster than either component system would on its own. 

We now give our explicit bounds on the distance to stationarity for the case B < B* . 

Theorem 7. For all B G {0,B*) and ai,02 E K, let a = max(B, \ai\, |a2|)- Then there exists 
]^B,ai,a2 s-t. n > NB,ax,a2 implies that for all t > 1, 

and 

IT.+a,v^,<rn+a.v^ W-^<rn+a.V^lMI < 4(i?M-^(l + n-|)eX^ 

Note that Theorem 7 provides a bound for any sufficiently large fixed n and all times greater than 
1, which is stronger than the bounds that would follow by naively applying the weak-convergence 
theory. Indeed, using the weak-convergence theory one could derive bounds of the form 'for any 
fixed time t, there exists a sufficiently large Nt s.t. n > Nt implies . . . ', but these bounds would be 
too weak to make statements about all time for any fixed n. We note that although we were able to 
derive partial results for the case B > B*, the derived bounds were considerably more complicated 
than those of Theorem 7, and we leave it as an open question to derive simple explicit bounds for 
the case B > B*. 
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We now discuss a practical 'rule-of-thumb' interpretation of Theorem 7 in the context of man- 
aging call centers. Suppose one is operating a call center in the H-W regime (with < B < B*) 
and at time t = the center begins either very overloaded (Q(0) = n + a^/n) or underloaded 
(Q(0) = n — a^/n) with a ^ a large constant (independent of n). Then Theorem 7 says that for 

any e > 0, the probability that a caller at time ^ = -^a"^ + ^log{^) has to wait for service will be 
roughly within e of the probability that a caller in steady-state has to wait for service. In particular, 
after t* ^ time units have elapsed, the system will have returned to within e of steady-state. 

3. Technical Preliminaries . In this section we introduce several notations, functions, and 
results from the literature that will be critical for our analysis. We then give a new characteriza- 
tion for 7„ (Lemma 3), which will be amenable to asymptotic analysis. We also prove some basic 
asymptotic results that will be used throughout the paper, and state some additional properties of 
the parabolic cylinder functions. 

We begin by associating several functions to the M/M/n queue, as in [29] and [38]. Let 

1 if k = 0; 

1 + \n — X if k = 1; 

{\n + k- x)fn,k~i{x) - \n{k - l)fn,k-2{x) Otherwise. 

By a simple induction argument, for < /c < n, 

(7) fn,k{^) = Y.{^)j^n\i{i-x). 

Also note that 



fn,k{x) = < 



(8) 



Z{fn,n-i) > 1, and /„,„_! (x) > for X < Z{fn,n-l] 



Indeed, for x < 1 and < j < n — 1, J|*L/(« — x) > 0, and by (7), fn,n-i{x) is a non-negative sum 
of such terms with the (n — l)st term {X^"^) strictly positive. Let 



a„(x) = i 



l(^Xn + n- X - ((A„ + n-xf - 4A„n)2^ if x < {y/n - ^/X^)'^ ; 



oo 



otherwise. 



Note that a„(x) is real-valued for x < {^/n — -v/A^)^5 aiid therefore 

1 



(9) 

Also, 
(10) 

We also define 



an{x) = ^{^Xn + n - X - {{^/n - ^/Xnf - x) 2 [{y/n + ^/Xuf -x)'^^. 



n. 



aoo(x) 



^{\[B-{B^-Ax)'2) ifx<f; 



oo 



otherwise. 
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Let 



bn{x) 



A j (4A„n- (A„ + n-x)2)2 if - ^/Kf < x < + ^/X 



2. 

nj ) 



I oo otherwise. 

Note that bn{x) is real-valued for x G [{\/n — (\/^+ VXi)^]) ^iid therefore 

^-l^-l^"! 1, t ^\ / /IT /\ \2\- ( f rz I A \2 



Let Cr„(y) = /„,n(y) - \/\^/n,n.~l(y), '^l^rXv) = /n,n(y) " an{y) fn,n-l{y) , 

( fn,k{y) 



Zn,k{y) 



if/„,,_i(y)/0; 
oo otherwise; 



and Zn{y) = Zn,n{y)- We also define 

a/^!3)+^ ifi?,-i(-5)^0; 



2;oo(y) 



oo 



otherwise; 



Let 



A Un.(y) - an(2/) if / OO or a„(?/) / oo; 



oo 



otherwise; 



and 



(12) Qn,k{x) = { 



An 



if k = 0; 
if k = 1; 



^ ~ ^ + ""° An '^'"'' )Qn,fc-i(a;) - ™°(^^ '^'"^ Qn,fc-2(a:) otherwise. 

We note that the Qn,k{x) polynomials can be associated with the well-studied Poisson-Charlier 
polynomials (see the discussion in [38]). It is shown in [38] that 

^ ^2 



(13) X e [{^/n - \/\nf, iVn + vKf) implies c„(x) = (5^_„(x) - Qn,n-l{x)Qn,n+l{x) > 0. 
Let 



9n{k) 



A 



( A„ \k—n 



if < /c < n; 
otherwise. 



We now cite several results from [29] and [38] that characterize 7„ and P^j{t). 



Theorem 8. [29] If Z (an) > (V^-^f, then = {V^-VKf- If Z{an) < (V^-^)^ 
then 7„ = Z(V'n). 

Theorem 9. [28], [38] If ^ > pi, then for all 
(14) Pr.it) - P;(oo) = (27r)-i5„(j) — (A„n)-i / e-"*Q„,,(x)Q„j(x)64x)c„(x)-idx. 



We now cite some properties of Zn{x),an{x), and ^ni^) stated in [29], all of which follow from 
elementary calculus and simple induction arguments, or the interlacing property of the zeros of 
orthogonal polynomials [8]. 

Lemma 1. [29] 

(i) Z{fn,k-i) > Z{fn,k) for k <n. 

(ii) For X G ( — oo, Z(/„,^„_i)) and k < n, Zn,k{x) is a strictly positive, continuous, and strictly 
decreasing function of x. 

(Hi) For X G ( — oo, (y^ — \/A^)^], an{x) is a strictly positive, continuous, strictly increasing 
function of x. 

(iv) For x G ^ — oo, min ((-^/n — -v/A^)^) -^(/ra.n-i))^ ; '^n{x) is cl continuous, strictly decreasing 
function of x. Also, if {y/n — \/\n)'^ < ^(/n,n-i); then ^'„(x) is continuous at {y/n — y/X^Y ■ 

We now prove some additional properties of {zn^k{x)}- 

Lemma 2. Zn^k{x) is a twice- differentiable concave function of x on (0, 1) for all k < n. 

Proof. We proceed by induction on k. For the base case, consider k = 1. By definitions, 
Zn^i{x) = A„ + 1 — X, and the assertion is trivial. Now, let us assume the statement is true for 
j = l,...,A:-lwithA;-l < n-l.By definitions, /„,fc(x) = (A„+/c-x)/„,fc_i(x)-A„(A:-l)/„,fc_2(x). 
Since by (8), Z{fn,n-i) > 1, we have by Lemma l.i that fn,k-i{x) > and fn,k-2{x) > 0, from 
which it follows that Zn^k-i{x) > and Zn,k{x) = (A„ + A; — x) — Xn^k — l){zn^k-i{x))~^ . Thus 

= -Xn{k - l)^^{Zn,k-l{x)~^) 

= -Xn{k - - Zn,k-l{xy'^-^Zn,k~l{x)) 

- 2Zn,k-l{x)~^ { — Zn,k-l{x)) + Zn,k-l{x)~'^ ^ Zn,k-l{x] 

Since Zn^k-iix) > 0, (^z„,fc_i(x))^ > 0, and by the induction hypothesis ^Zn,k-i{x) < 0, it 

follows from (15) that z„^fc(x) is twice-differentiable on (0, 1) and satisfies ^^Zn^k{x) < (concavity), 
proving the induction. □ 

We now prove a slight variant of Theorem 8, which will be more amenable to asymptotic analysis. 

Lemma 3. 



(i) If {y/n - y/Xi)"^ < 1 and ^n{{Vn - y/Xnf) < 0, then Z(^'„) = Z+(^'„) is the unique zero of 

\I'„(x) in the interval (O, {y/n — -v/A^)^) and 7„ = Z+(^'„). 
(ii) If{y/^-VKf < 1 and ^„((V^ - \/fc)^) > 0, then = {y/n - VKf ■ 
(Hi) If {y/n — \/A^)^ > I, then Z{"in) = ■2^^(^n) is the unique zero of'^n{x) in the interval (O, l) 

and 7„ = Z+(^'„). 
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Proof. We begin by proving some properties of and (T„. From (7), 



V" ('^\ T-rn-fc . 
2^fc=0 fc! 



1 



- A„, + n-((A„ + n) - 4A„n) 



> 0. 



If i^ffi-^,)^ > 1, then 



*n(l) 



ELoa)Amr=t(^-l) 1 



Xn + n-1- ((A„ + n - 1)2 - 4A„n) 



An + n - 1 - ((A„ +n-lf - AXnu) 



= -(^A„-n + l+ ((A„-n + l)2-4An) 

< ^(A„-n + l + |A„-n + l|) 
= since A„ — n + 1 < 0. 



We have 



k=0 



[n 



n— 1 / n—l—k 

n — 1 



1=1 



k=0 



n 



1=1 



k=0 



k=0 

A;! 



k=0 



> 0, 



and 



n—k n—1 / ^\ n—l—k 

-nil) = Ed^-iK^-D-v^ErrV" n(^-i) 



fc=0 



fc=0 



4 = 1 



A" — \/ XnnX„ 



n-l 



< 0. 



We first prove assertion i. Since {^/n—^/X^^)'^ < 1, we have by (8) that min i{^/n—^/X^)'^ , = 
{^/n—^/X^)'^ , and fn,n-i{x) > on ( — oo, {\/n—y/\l)'^]. Thus from definitions (and dividing through 



by fn,n-i{y) ), we find that is the same sign as Zn{y) — V^nn, and that ipn{y) is the same 

sign as ^'n(y), on ( — oo, {^/rl — \/A^)2]. Also, by Lemma l.ii, Zn{y) is strictly positive, continu- 
ous, and strictly decreasing on ( — oo, (-^/n — \/A^)^]- Thus by the Intermediate Value Theorem, 
ajiiy) has a zero on ( — oo, {^/n — -v/A^)^) iff •2n((\/^ — — VKin < 0. By definitions and 

(10), we thus have that Zia-n) < {Vn — V%i)^ iff ^n((v^ ~ VXi)^) < 0. Since by assumption 
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^n((\/^~ < 0, it will be the case that Z(o"„) < {\/n — -v/A^)^- Furthermore, it follows from 

the fact that ^'n(O) > and the continuity and monotonicity of ^n{y) guaranteed by Lemma l.iv 
that Z{'^n) will be the unique zero of ^'n in the interval (O, {\/n — \/A^)^)- Since ipniu) is the same 
sign as ^'n(y) on ( — oo, (v^ — VKi)'^], we will also have that Zi^^) = ^(^n)- This, combined with 
Theorem 8, completes the proof of assertion i. 

We now prove assertion ii. Since {y/n—y/X^)'^ < 1, we have by (8) that min {{y/n— y/Xnf' , -^(/n.n-i)) 
{^/n— y/'Kif' , and fn,n-i{x) > on (— oo, (-v/n— \/A^)^]- Thus from definitions (and dividing through 
by fn,n-i{y) )) we find that an{y) is the same sign as Zn{y) — \/ A^n, and that ipn{y) is the same sign 
as ^'n(y), on ( — oo, (-^/n — -v/A^)^]- Also, by Lemma l.ii, Zn{y) is strictly positive, continuous, and 
strictly decreasing on ( — oo, {y/n — \/A^)^]. Thus by the Intermediate Value Theorem, an{y) has a 
zero on ( — oo, (v^~ V%i)^) iff •2„((\/n — \/A^)^) ~ VKj^ < 0. By definitions and (10), we thus have 
that Z{an) < {y/n- \/A^)^ iff ^n{{Vn- y/K.)"^) < 0. Since by assumption ^'„((^- \/A^)^) > 0, 
it will be the case that Z{an) > (\/^~ V^)^- This, combined with Theorem 8, completes the proof 
of assertion ii. 

We now prove assertion iii. By assumption, (-^/n — y/X^)'^ > 1. Thus since Z(/„^„„i) > 1 by (8), 
we have that fn,n-i{x) > on (— oo, 1]. Thus from definitions (and dividing through by fn,n~i{y))-, 
we find that ipniy) is the same sign as ^'n(y) on (— oo,l]. Since cr„(y) is a polynomial and thus 
continuous, by the Intermediate Value Theorem and the fact that cr„(0) > and cr„(l) < 0, we have 
Z{an) < 1- Thus Z{an) < {^/n — By the continuity and monotonicity of ^'n(y) guaranteed 

by Lemma l.iv and the fact that ^n{^) > 0, ^'n(l) < 0, we will have that Z{^n) is the unique 
zero of in the interval (0, 1). Furthermore, since ipniy) is the same sign as ^niy) on (— oo, 1], we 
will have that Z{ipn) = Z[^n)- This, combined with Theorem 8, completing the proof of assertion 
iii. □ 

We now prove some simple asymptotic bounds for use in later arguments. 



Lemma 4. -| < ^/n — -y/A^ — "f "I ^ — ; '^"■^ hmn-s.oo(v^ ~ VKi) 



4n^-2B 



B 
2 • 



Proof. 



n 



(16) 

Note that 



B 



l + (lf)^ 



B 



l + (l-Bn-h¥ 



1 - Brr2 < (1 - Bn'2)2 < 1. 



Combining with (16) we find that y/n — -v/\i is at least and at most 



(17) 



B 



B 



2 

B^ 



2 - Bn-2 

The lemma follows letting n — )■ oo in (17). 



4n2 - 2B 



□ 



A _ 1 

Let qn = 1 + n 4 . 
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Lemma 5. For each B > 0, there exists Nb s.t. n > Nb implies gn{k)^ < ^n^^ for all k 



Proof. Since ^ < 1, we have that gn\k) < j^— |T foi' all A; > 0. Thus 
n! Ak^"J n! 



5n(A;)— < 



n2 — i? 

< exp(^ + — ) 

n2 — B n2 — B 

= exp{B^ + — + — ) 

n2 — B 712 — B 



e-^%xp (0(n 2)) 



< Qn^^ for all sufficiently large n by a simple Taylor series expansion. 



Lemma 6. For all < x < I, T > 3, n > T + 1, and k e [0,n - T - 1], 

n—k 

(18) (n - A;)-^r^ < H {1 - -) < e^^~' {n - k)-^T\ 

j=T+l ^ 

Proof. Since < x < 1, by a simple Taylor series expansion we have that for j > 3, 1 < 
1 + jz. Thus for j>T + l, 

n—k —2. n—k 



n s n (i+i 

j=T+l ^ j ' j=T+l •' 



1 

1^ 



< n 

j=T+l 

EOO 1 



and 



n—k n—k n—k 

e J . 



(19) e-- n ^-^^ n (1-7)^ n 

j=T+l i=T+l j=T+l 

From [45], we have that for all k > I, log(fc) + C < ^^fc < ^ + log(fc) + C, and thus log(^) - ^ 
Hn-k -Ht< log(^) + 

2(^n-k) ■ Furthermore, since by assumption n — k > T, we have 

/ X , ,n — k. 1 , ,n — k . 1 

(20) log(^) - ^ < Hn^u -Ht< log(^) + — . 
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Combining (19) and (20) yields 



n—k 

(21) e-^"'e-#(n- A:)-^r^ < J| (1 - -) < e# (n - A;)-^r^. 

i=T+i ^ 

Since < x < 1 we have < ^ < and the lemma thus follows from (21). □ 

We now state some additional properties of the parabolic cylinder functions. Dx{z) is an entire 
function of z for fixed x ([15], Chapter 8, Section 8.2) and an entire function of x for fixed z [9]. It 
is stated in [21] that for all G M, 

(22) D^+i{z) - zD^iz) + xD^_i{z) = 0; 

(23) ^^D^{z) + ]^zD^{z)-xD^-i{z) = d. 

4. Asymptotics of 'J'^(a;) . In this section we derive the asymptotics of '^n{x) and ((\Ai~ 
•\/A^)^). In particular, we prove 

Theorem 10. For all B > and x G (0,1), x < lim„^oo = ^00(2;). 

We also prove 

_ 1 

Corollary 2. For < B <2, lim„^oo A„ ^^'„((^- = C{B). 

We proceed by separately analyzing the asymptotics of and fiiif^^^ and then use the 

fact that by definition - £iiM^ _ M£)^ 



We now analyze the asymptotics of an{x 

B'^ 1- an(x)-X 



(24) 



Lemma 7. ForO<x<^, lim^^oo = aoo(x). 

^ V An 

Proof. Prom (9) and Lemma 4, 

an{x) - \n _ By/n - X - {{y/u + \/A^)^ - x) ^ ((^ - ^/X^)'^ - x) 



2x/a: 



The lemma then follows from the fact that lim„_>.oo '^J! ^ = lim„_i.oo ^^^'^'^f^ ^ — = 1, and 

•2V An ^ .2V An 

Lemma 4. □ 
We now analyze the asymptotics of Zn (x) , proving that 

Proposition 2. For x g (o, l), lim„^oo "'"^"j^^" = .200(2;). 
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By (8), fn,n-i{x) > on (0, 1), and thus 



(25) 



V Ar), 



VA;E^=o(V)A^n-=i'-'^o---) 



The remainder of the proof of Proposition 2 proceeds by a series of lemmas. 
Lemma 8. 

. , . I"-^n - 51 p-^"^" 

(26) 



Znix) -Xn_ /^ELo(- - ^) n (1 - f) e-^" t 



^A„ 



ELo^n-=f(i 



Proof. We begin by analyzing the numerator of (25), which is equivalent to 



n " / \ n—k n— 1 

.7 = 1 k=l ^ ^ 1 = 1 



fc=0 



n — 1 
(fc + 1)- 



n-{A;+l) 



^)a^^ n (i-x) 



(27) 



n 

no 

n 

no 

n 

no 

n 

no 



fe=i 

n 



n—k 



fc=l ^ ^ 



'n- 1 



k=l 
n—k 

[])>ij{{j-x) 

+E|!^("-'M^no-) 

fc=i ^ ^ i=i 



n—k 



+ (n- l)!^(n- A;) 



n-=f(j-^) Aj 



A;=l 
n 



+ (n-l)!5^(n-A;)n(l-7) ^ 

fc=l jr = l 



(n-A:)! A;! 

n—k 

^ ~ j' k\ 



n—k 



{n-l)\{nJ{{l--) + Y,{n-k)\[{l 



J 



n—k 



k=l 



X A: 
7^ ^! 



n 



i)!E(--^)n(i-7)| 

fc=o i=i 



We now analyze the denominator of (25), which is equivalent to 



n-l 



A;=0 



n — 1 
(fc + l)-l 



ifc+i 



n-(fc+l) 

n 



n , / \ n—k 



k=l 



n n—k 



(28) 



A;nn-l)!E^n(l-7) I- 

fc=o j=i ^ 
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Combining (27) and (28), and normalizing both by ^^-1)1 ' yields the lemma. □ 

Note that by a simple application of Stirling's Inequality, we find that lim„__>oo Vne " "[a~|T ~ 
{2^/^^)~^. Thus we may fix a positive integer T > 3, and assume that n is sufficiently large to ensure 
that ^ < 2, [n- r-^V^ + 1 < n - (r+ 1), ^ < 57, and e"^"!^ < n's. 

Lemma 9. The numerator of (26) (ignoring the \/A^ prefactor) is: 
(i) hounded frova below by 

T n-(T+l) n-k 

(29) Yin-^) E (--^) n (1-7)^'' F' 

j=l A:=0 j=T+l •' 

(a) hounded from above by 

n-(T+l) n-k 



(30) 



j=i ^ k=o j=T+i ^ llj=i(i - yJv?^ 



Proof. Note that the numerator of (26) is equivalent to 
(31) 

T ,n-{T+l) n-k \k T n n-k , ^ x 

n(i-£)( Y. n e (»-')n(i-7)=-'"|)' 

j = l ^ fc=0 ^=^+1 ■ i=l k=n-T j=l ■' ' ^ 

The second summand in (31) is non-negative, yielding a trivial lower bound. We now obtain an 
upper bound on the same expression. Since the mode of a Poisson r.v. with mean is [A„J , and 

by assumption e " "[a~|T — ^ ^ , we have that for all k > 0, 

\k \l^n\ 

,32) e-.|<,-..^<„-. 

We now use the fact that n — k < T + 1 for k > n — T, 1 — j ^ Ij and (32) to conclude that the 
second summand in (31) is at most 

(r + i)2 



(33) 



n?=i(i-f)^/^' 

Combining our lower and upper bounds on (31) completes the proof of the lemma. □ 

Lemma 10. The denominator of (26) is: 
(i) hounded from below by 

T n-(T+l) n-k 

(34) n(i-7) E ^ n (i-7)^"'"F' 

j=l k=0 j=T+l •' 



15 



(ii) bounded from above by 

3=1 ^ k=0 3=T+1 ^ llj=l(l-j) 

Proof. Note that the denominator of (26) is equivalent to 



T Ji-(T+1) n-k T n n-k 

(36) n(i-7)( E ^ n a-7)^-^4+(n(i-7))" e ^na-7)^ 

j=l fc=0 j=T+l •' ' j=l •' k=n-T j=l 



We now bound the second summand in (36). Note that it is non-negative, yielding a trivial lower 
bound. Also, it is at most 

max A:|T(1-- e-^"-y < n— by 32 

nJ=i(i-f)"-^^'=^" M ^ ^' " nj=i(i-f) ' ' 

(37) < i^±ll^. 

" nj=i(i-f) 

Combining our lower and upper bounds on (36) completes the proof of the lemma. □ 
LetCT,. = max(^A^^,^^^_^). 
Lemma 11. 



(38) 



II 



(39) 



_ 1 

r.n 2 



At- 



-2T-1 v^n-CT'+l) k ( 



^ kl 



Proof. Combining Lemma 9 and Lemma 10 we obtain 



(40) 



Xr 



> \ X 



k\ 



Sfc=0 



'(^-fc)n"=T+i(i-f) 



^n-(T+l) 



fcn-=T+i(i 



j) ^ fc! + 



n=i(i-f) 



and 



(41) 



Zn{x) - Xr 

\/a7 



E 



n-(T+l) 
A:=0 



(^-^)n-=T+i(i-f)e 



k\ 



+ 



(r+i)2 



< VA, 



n;=i(i-f)v^ 



z^fe=o '*'lli=r+ivJ- , J e k\ 
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Combining (40) and (41) with Lemma 6 and the assumptions on n, we find that 



Xr, 



k\ 



1 „_(T+1)_^ _^ (T+l)T-^ ' 

^ ^k=o v^l" f') e fc! + V^"n-=i(i-f) 



from which (38) follows since < 2. (39) follows from a nearly identical argument, proving the 



lemma. 

Lemma 12. 



□ 



(42) 

■2-n(3^) Xn 



> 



k\ 



^ ^n-T-^V^ k{n-k)-^ e-AnAfe ^ 2e^^-^ T-d-^) ( j^) ^ £^ 



fe=0 



+ 



1-x 



(43) 



-2T-1 Y^rn-r-l0i] k{n-ky e-^"Afc 



Proof. We first prove i. The term ^Y^=q"~^^^ ^ ^ ^"iJ appearing in the denominator 

of (38) is at most 



E 



fc=0 



(44) 



< 



< 



k{n - k) ^ e '^"A„ 
~k\ 



vAn 

s E 

A:=0 

s E 

fe=0 

s E 

A:=0 

E 

A:=0 

E 

fe=0 



+ max 



\n-T-^^/Ti\+l<k<n a/A. 

k{n - ky^ e'^^Xn n 



A;=[n-T-iv^ + l 





/ An 




k\ 


fc(n 


-fc)- 


e- 


-An \ fc 




/a; 




/fc! 


k[n 


-k)' 


g- 


-An \fc 








k\ 


k[n 


-ky 


e- 


-An \fc 








/fc! 


k{n 


-k\ 


g- 


-An \fc 




J An 




/fc! 



H ^=''^ 2 

An 



+ 




Ar 



^ (n-A;)-- by (32) 

fc=|"n-T-iv^+l 
[T-ivM-l 

j=r+i 



+ 




+ 2 — 

An 1-x 




n JO 



n 



n 



H since — < 2. 

1-x Xn 
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We may thus upper-bound the denominator of (38) using (44), and then multiply both numerator 
and denominator by An^ , to prove i. 

We now prove ii. By non-negativity the denominator of (39) is at least e~^'^ ^ Y^k=^ ^^"^^'^ — ~~H~^' 

3;-! 

and thus ii. follows from (39) and multiplying both numerator and denominator by A^^ • □ 

We now relate the summations appearing in (42) and (43) to the expectations of certain functions 
of a poisson r.v., and then show that these expectations converge to certain integrals as n — )• oo. 
Let Xn denote a Poisson r.v. with mean A„, Z„ = ^ and I{G) denote the indicator of the 

event G. Thus 



(45) Xn — V A„Z„ + A 
and 

n- Xn = n- {■s/Xn'Zn + A„) 

(46) = B^-^nZn. 
Furthermore, let 

(47) n,n = (i?yj - Znf-AZn < B - ^), 

(48) Y,,n ^ Zn{B^ll-Zn)-^IiZn < (B - T-^) ,/f + ' li" " ) , 

V '^n, V '^n V An 

and 

(49) Y,,n ^ (B^^ - ZnT^HZn < {B - T-)^^ + \-T-'^^n-T-^V^) ^_ 
We then have that 

Lemma 13. If \min^aoE\Yi.n\, lim„_>oo -E'[l2,n]; and Ivain^oo E\Y'i^n] M exist, and are finite 
and strictly positive, then 



) — A„ _4T-l 1™ 

(50) limmf -= > e 



hm^^oo E[Y:,^n\ + ^^^^^ ' 
(ii) 

(51) hmsup < e^T- lim„^oo E[Y,,n] 



/A^ lim„^oo-E'[>3,nJ 



Proof. We begin by expressing the summation appearing in the numerator of (42) and (43) as 
an expectation. 



\y Yl (n-fc)^-"e-^"^ = Xn''E[{n-Xn)^--l{Xn<n-{T + l))] 



k\ 

k=0 



= Xn' E[{BV^-^/X;,Zn)'-^I{^/KZn + Xn<n-iT + l))] by (46) 
(52) = E[Yi,n]. 
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We now express the summation appearing in the denominator of (42) and (43) as the sum of two 
expectations. 

= A^^[X„(n - X„)-"/(X„ < [n - T-^V^l)] 
by (45) and (46) 

(53) = \7^E[{Zn + ^/K){BV^-^/KZn)-''I{y/KZn + \n<\n-T-^V^)]. 



We now manipulate the indicator variable appearing in (53) to make it more amenable to asymptotic 
analysis. The indicator is 



I{Zn < 
I{Zn< 



\n — T~^^/n] — (n — B^/n) ^ 

\n-T-^^ -{n-T-^^) + {B -T-^)^). 

V An, 



(54) = /(Z„ < [B - T-\[^ + - li" - ). 

V An V An 

Plugging (54) back into (53), we find that the summation appearing in the numerator of (42) and 
(43) is 

An-ii;[(Zn + V^)(i^,/f - Zn)--I(Z„ < {B - T-)./f + " ^^^^ li" " ^^^^ )] 
V '^n \ vAn 

(55) 

_ 1 

= An ^ i?[l2,n]+-S[y3^n] by breaking up the (Z„+y An) term and applying linearity of expectations. 

The lemma then follows by plugging (52) and (55) into (42) and (43), and observing that 

lim = , lim A„2 Ctx = 0, and lim A„2^[y2n] = 0. 

n— >oo 1 — X 1 — X n^oo ' n— >cxd ' 

□ 

We now explicitly compute the limiting values of all expressions appearing in (50) and (51), allowing 
us to compute lim„_s.oo . 

Lemma 14. 



(56) 



n^oo VAT rB "'^ 



j^^{B-zY-e--dz 
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Proof. It is easily verified that {Yi^n},{y2,n},{Y3,n} are all sequences of uniformly integrable 

r.v. Let /i(y) ^ {B - yy-^I{y < B),h{y) = y{B - y)--I{y <B- T~^), ^y) = {B - y)-I{y < 
B — T~^). Let denote a normal r.v. with mean and variance 1. It may be easily verified that 
{^i,n}, {^2,n}, {^3,n} Converge weakly to fi{N), f2{N), /^{N) respectively, where all three conver- 
gences follow from the well-known convergence of the scaled Poisson r.v. Zn to the Normal r.v. N. 
It thus follows from the uniform integrability of all three sequences that 

(57) lim E[Yi,n] = E[h{N)] = {27t)-2 / {B - zf^'^e-^dz, 



(58) lim E[Y2,n] = E[f2{N)] = (2^)-2 / z{B - z)-'^e--dz, 

-/ — oo 

and 

(59) lim ElY^^n] = E[h{N)] = (27r)-2 / {B - z)~^e-^dz. 

-/ — oo 

The lemma follows from plugging the limits (57)-(59) into Lemma 13, realizing that T was general, 
and letting T — )■ cxd. □ 

We now relate the integrals appearing in (56) to the parabolic cylinder functions. In particular, we 
prove 



.2 



Lemma 15. ^ — = Jr , L + B. 

Proof. Note that the right-hand side of (56) is equivalent to 



t-e-'-^dt e-^ t-e-'^+^'dt 

^oo^_((._2)+l)^_4_(_B)t^^ 
^oo^-((._l)+l)^_^+Bt^^ 

f,^t-i^--^^+')e-'^-^-^^'dt 

r{-{x-2))D^_2{-B) . 
= —7 by (2), smce x G (0, 1) 

_ .{D^hB)+BD^-i^B)l 

(60) = ^^'*'7^j>x +B since r(2- x) = (1 -x)r(l -x). 

E>x-i{-B) 

□ 
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We now complete the proofs of Proposition 2, Theorem 10, and Corollary 2. 

Proof of Proposition 2. Proposition 2 follows by combining Lemma 14 and Lemma 15. □ 

Proof of Theorem 10. Since by definition = — a„(j;), Theorem 10 follows im- 

mediately from Lemma 7 and Proposition 2. □ 

Proof of Corollary 2. Suppose < e < min(^, 1 — -^). Then for ah sufficiently large n, 
< — e < (\/n~\/\i)^ < min(l — e, ^ + e) by Lemma 4, and < 1 — e. By Lemma l.ii and (8), 
we have that ^"^"^l ^" is a continuous, strictly decreasing function of x on (— oo, {\/n— -v/A^)^ + e]- 



It follows that for all sufficiently large n, ^"^ < "(^^ ) — " < ^"^ ^ ^" _ Thus by 

Proposition 2, for all sufficiently small e > 0, 

(61) 

Zoo{—r + < limmt = — < hmsup = — < Zoo{— e). 

4 n->oo ^/Xn n^oo V K 4 

To proceed, we will now prove that Zoo{x) has certain continuity properties. In particular, 

(62) ^00(2;) is a differentiable, non-increasing function of x on (—00, 1]. 

Indeed, by (4), Dx-i{—B) > for x < 1. The differentiability of ^00(2;) on (— 00, 1] then follows 
from the fact that Dx{—B) is an entire function of x. That Zoo{x) is non-increasing follows from 
Lemma l.ii, (8), and Proposition 2. 

The continuity of 2^00(2;) in a neighborhood of (guaranteed by (62) since < ^ < 1) , along 

with (61), thus implies that lim„_j.oo — 7y= — = •2oo(-x)- (^0)' — ^ ~ 



— VKi- Thus by Lemma 4, lim„^oo ^/^^ ^ — ~ = "f ; ^^'^ the corollary follows since 

C(i?) = ^oo(f ). " □ 

5. Asymptotics of Z~^{'i/n). In this section we derive the asymptotics of Z+(^'„). In partic- 
ular, we prove 

Theorem 11. If B <2 and C{B) < 0, or B>2, then lim„^oo 2'+(^'„) = Z+(^oo)- 

We begin by establishing some additional properties of ^oo{x)- 

Lemma 16. ^00(2;) is a continuous function of x on ( — 00, min(l, ■^)) , left- continuous at 
min(l, ^), and a differentiable function of x on ( — c«,min(l, -^)), satisfying -^^ooix) < —{B'^ — 
4x)~^ /or X < min(l, ^). 

Proof. That aoo{x) is differentiable on (—00, ^), left-continuous at and satisfies ^000(2;) = 
(B^ — 4x)~^ for X < follows from elementary calculus. The lemma then follows from (62) and 
the fact that ^^00(2:^) = Zoc{x) — aoo{x). □ 

We now prov6 some ctdditioncil properties of ^"'"(^^o). 
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Lemma 17. 



(i) If B <2 and C{B) < 0, or B>2, 

(a) ^oo{^) has a unique zero Z^{^ao) ^ (O, min(l, ■^)), 
(h) ^-ooW >0 on [0,Z+(^oo)), 
(c) ^'oo(x) < on (Z+(^'oo),mm(l, ^)]. 
("ii; If B <2 and ({B) = 0, 

(a) ^ooix) > on [0,min(l,^)), 

(b) Z+{^^) = f. 

Proof. We first prove assertion i. By Lemma 16, ^'00(2;) is a strictly decreasing function of x 
on [0, min(l, -^)) , continuous on [0, min(l, -^)) , and left-continuous at min(l, ■^). ^'oo(O) > 0, 
since ^^{0) = ^2i^ + B, and by (3), ^2l_^ > 0. 

Also, ^'oo(min(l, ^)) < 0. Indeed, if B > 2, min(l, -^) = 1. Furthermore, ^'oo(l) < 0, since 
^oo(l) = + 1{B + {B^- 4)-2), and by (3), 

^i|^ + l(i? + (B2_4)l) = + 1(5 + (5^-4)1) 

= i((i?2_4)i_5) 
< 0. 

If S < 2, then min(l, = and ^'oo(^) = C{B) < 0. Combining the cases 5 < 2 and B > 2 
demonstrates that ^'oo(™in(l, ■^)) < 0. 

Thus ^' 00(0) > 0, and^'oo(min(l, -^)) < 0. That 00 (x) has a unique zero 00) G (O, min(l, ^)) , 

^ooix) > on [0,Z+{^^)), and ^'oo(x) < on {Z+{^^) , min(l, ^)] then follows from the conti- 
nuity and strict monotonicity of ^ao{x) on ( — oo,min(l, ■^)] guaranteed by Lemma 16. 

We now prove assertion ii. The lemma follows from the fact that ^'00(2;) is strictly decreasing on 
( — oo,min(l, ^)) and left-continuous at ^ by Lemma 16, and C{B) = 0. □ 

Proof of Theorem 11. We first show that for all sufficiently large n, is the unique 

zero of in the interval (0, min(l, {\/n — V%i)^)- Indeed, if S < 2 and C{B) < 0, then for all 
sufficiently large n, min(l, (\/n — -v/A^)^)) = {V^~ V%i)^ by Lemma 4 and ^'„((y^— -v/\i)^) < 
by Corollary 2. It follows from Lemma 3.i that for all sufficiently large n, is the unique 

zero of in the interval min(l, {^/n — If instead B > 2, then for all sufficiently large n, 

{^/n — a/A^)^ > 1 and min(l, (-v/n~ = 1 by Lemma 4. It follows from Lemma 3.iii that for 

all sufficiently large n, Z+(^'„) is the unique zero of in the interval (0,min(l, {^/n — -v/A^)^))- 

Now suppose for contradiction that liminf„_>oo Z~^(^n) < ■^"'^(^oo)- By Lemma 17. i, Z~^(^oo) £ 
(0, min(l, -^)) , and ^'oo(a^) > on [0, Z+(^'oo)) . It follows that there exists e > such that 
< liminf„_^oo -^^(^n) + e < min(l, ^), and ^oo(liminf 

n^oo Z^{'^n) ~l" > 0. Thus by Theorem 
10, for all sufficiently large n, ^'n(liminf„_j.oo ■Z^^(^n) + e) > 0. By the monotonicity of ^n{x) (from 
Lemma l.iv), it follows that for all sufficiently large n, ^'n(a^) > on (— 00, liminf„_^.oo Z^{'^n) + (-)- 
But by the definition of liminf, there exists an infinite strictly increasing sequence of integers {ui} 
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s.t. Z+(^'„J < liminf„_^oo ^"^(^n) + e for all i. Thus for all sufficiently large i, ^'„^(Z+(^'„.)) > 0. 
This is a contradiction, since from definitions = for all i. 

Alternatively, suppose for contradiction that limsup„_^oo > Z+(^oo)- By Lemma 17. i, 

Z^{^oo) e (0,min(l, ^)), and ^'oo(a:) < on (Z+(^oo), min(l, ^)]. It follows that there exists 
e > such that < limsup„__^o^ — e < min(l, -^), and ^'oo(linisup^_i.oc — e) < 

0. Thus by Theorem 10, for all sufficiently large n, ^ri(limsup„_j>oo ■^^(^n) — e) < 0. By the 
monotonicity of ^n{x) (from Lemma l.iv), it follows that for all sufficiently large n, ^n{x) < on 
(limsup^^^ Z+{^n) - e,min(l, {y/n - y/Xnf)) . 

By the definition of limsup, there exists an infinite strictly increasing sequence of integers {rii} 
s.t. Z^[^n.) > limsup„_^o^ — e for all i. Combining the above, we find that for all sufficiently 

large i, < 0. But this is a contradiction, since from definitions = 

for all i. We conclude that for the case B < 2 and C(-S) < 0, or i? > 2, liminf„_j.oo = 
limsup„^^ = Z+(^oo). □ 

6. The Zeros of C(-^)' In this section we characterize the set of € (0, 2) for which C(^) ^ 0, 
which will allow us to use Theorem 11 effectively. In particular, we prove 

Theorem 12. B* e [\/2,2). ({B) > on [0,B*), and ({B) < on {B*,2]. 

We also complete the proof of Proposition 1. We begin by proving some properties of the derivatives 
of -^^^^5 which will then enable us to analyze -^C,{B) using the multivariate chain rule. 

Lemma 18. 



(i) Dx-i{y) G (0,oo) for X < 1 and ally. 

'■vL 
i{y) 



(ii) j^"" ^^)-. is a smooth function of y for x < 1, and 



(63) ^_54L., + (_54L)=_,.«^fe) 



dyD^^i{y) -Dx-i(y) 
(Hi) j!^"" ^^^y^ is a smooth function of x for x <1 and all y. 

Proof. Note that if x < 1 then Dx-i{y) > for all y by (4), from which assertion i. follows. 
We now prove assertion ii. Since Dx{y) is an entire function of y for each fixed x, it follows from i. 



that for X < 1, „ / n is a smooth function of y, and 



ID Ou OLLLVjyjhLL ±Ui±UtiWl- 

d D,{y) 



dyDx-i{y) 
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D.-i{y)iyD.{y) - D.{y)iyD.-i{y) 

D^.i{y){xD^.i{y) - \yD^{y)) - - l)D^-2{y) - \yD^-i{y)) 

xDl-M - \yDAy)Dx^i{y) - ~ i)D^{y)D^^2{y) + \yaMDx-i{y) 

DU(y) 

{x - l)D,{y)D,^2{y) 



by (23) 



X 



X + 



DLiiy) 

{X - l)D,{y){£^M^i£^] 

DLiiy) 

D,iy){D,{y) - yD,^i{y)) 



"2 by (22) 



Dl-iiy) 



x + { D^(y)f_y_D^iy) 



'^.-i(y)' "D^^iiy)- 
Also, assertion iii. follows from the fact that Dx{y) is an entire function of x for all y, and (4). □ 

Db2{-B) 

We now prove some properties of -dBlT~^~ipB) applyiiig multivariate chain rule to Lemma 
18. 



Lemma 19. 



(i) Dg^_^{-B) G (0,00) on (-00, 2]. 

r**/' D 'T' (- 'm ^ differentiahle function of B on (— cxd,2], and on (—00, 2] we have 



'b2 S-B) 

J- 



dBDg^ (-B) 2dxDx.i{y) 



.b2 



/ 4 \l 



D^{-B) 
B- * 



4 ^Db2 {-b)' {-by 

y=-B T ^ -4, i 



Proof. D ^2 ^{—B) > for i? < 2 by (4). Thus we may apply the multivariate chain rule to 

dJ(-b) 



find that ^ ~^ , — ^ is a differentiable function of B on (—00, 2], and 



d 



D^i-B) 



dBD, 



-B) 
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d B\ d D^iy) 



B d D,{y) 



+ 



^-^ dB 4 ' dyD,^i{y) 

y=-B 

d D^{y) 



-^iB^-^^ 



y=-B 



2 dxD^^iiy) 
B d D,{y) 



dyD,.,{y) 

y=-B 



y=-B 



2 dxD^.^i{y) 



y=-B 



„2 Ds2{-B) Ds2{-B) 
4 ^Z)s2 A-B)' (-B) 



by Lemma 18. ii. 



□ 



Although it seems clear from Figure 1 that C(-S) is strictly negative on {B*,2), the formal proof 

. Dg2{-B) 

-. Due to 



of this fact is somewhat involved, and will necessitate careful arguing about j^jy 

the fact that ^Dx{y) has no simple analytic form ( as opposed to -^Dx{y)i see (23) ), to proceed 

we will have to derive good bounds for ^ o'^^^^y) ■ bounds will rely on the concavity of certain 
functions. We begin by proving that 

Lemma 20. For any fixed y < 0, j^'^^^^y^ is a concave function of x on (0, 1). 

Proof. For any fixed B > 0, Zoo{x) is a concave function of x on (0,1), which follows from 
Lemma 2, Proposition 2, and the fact that pointwise limits of concave functions are concave. The 
lemma follows, since j^'^^j^^b) ~ -^oo(-^) ~ ^ O^y Proposition 2). □ 

We now use the concavity of ^^^^-^ and the associated monotonicity of ^7^^^^ to prove that 



Lemma 21. For any y < and xq € (0, 1), 

d DM 



dx Dx-i{y) 



< 



1 D^oiy) 



x=xo xoDxQ-i{y) 

Proof. By Lemma IS.iii and the Mean Value Theorem, there exists CxQ,y G (0,xo) such that 



d Dxiy) 



Thus 



(64) 



d Dx{y) 



dx Dx-iiv) 



< 



dx Dx-i{y) 
d Dxiy) 



x=xo 



dx Dx-i{y) 



D^oiy) _ Dojy) 



by Lemma 20 and the def. of concavity 



< 



Xo 

Xo Dxo-i{y) 



since > by (3), and xq > 0. 

D-i{y) 



□ 
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We now plug our bounds on -^jy^^^^ from Lemma 21 into our bounds on ^ ~^ from 
Lemma 19 to prove that 

Lemma 22. For < S < 2, 

, D^2{-B) 2 D^2{-B) „2 Ds2{-B) . 



dBD^i (-B) ^B 'Db2_ (-B) 4 ^Dg2^ (-B)' - B^ 

i ^ 4 ^ 4 ^ 

Proof. By Lemma 19. ii and Lemma 21, for < S < 2, 

, Ds2{-B) „ . D^2{-B) „2 Ds2{-B) D^2i-B) 

" — , ^ ~ ^ ~ \2 _ TJ ~ 



dBD^_(-B) '^^D^(-B) 4 'D^_(-B)^ D^_(-B) 



4 " 4 



^B 'd^ (-B) 4 "-D^A-BY ■ 

4 4, ^ 

This proves the first inequaUty. Furthermore, it follows that there exists G M (we may take 

^s2(-B) D^^(-B) 
^B = D ^ (-- Bj-^ such that jBTT^r~(FB) — ("B ~ B)xb — Xb ^- Since by elementary calculus 

the function (of x) (-^ — — — ^ attains its maximum at x = -^-^ — , we find that for all 
BG (0,2), 



dBD„2 A-B) - ^B 2 ' ^ 2 ' 4 



□ 



We are now in a position to study the zeros of C,{B\ 
Lemma 23. C(^) > on [0,5*), and G (0,2). 
Proof. By (3), 

and 

w X ^i(-2) -2e-i 
(66) C(2) = ^ + 1 = _ + 1 = -1<0. 

By Lemma 19. ii and definitions, C(-S) is a differentiable (and continuous) function of B on (— oo, 2]. 
It follows from the Intermediate Value Theorem that B* G (0, 2), and C,(B) > on [0, B*). □ 

Lemma 24. B* G [\/2,2). 
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Proof. It is proven in [40] (Theorem 4.1, i.) that 
(67) 7n > inf ( A„ + min(A;, n) — y/X^(^\/mm{k — 1, n) + \/ min(A;, 



k>l 



n, 



Note that for 1 < A; < n, 
An + min(A;, n) — \/A^(-\/ min(A; — 1, n) + min(A;, n)) 



A„ + A; - ^/Ki{y/k - 1 + \/A^) 
"v/A^ 



> - 



(68) 



> - 



2 ^/k 
Ia/AT 



k? + 



^/k^^fk^ 



For all /c > n + 1, the right-hand side of (67) equals {^/n 
and (68), we find that 



2 y/n 

A„)^. Combining this fact with (67) 



7„ > min(^^^^,(Vn- \f\'nf). 



(69) 



By Lemma 4, lim„_!.oo(\/^ ~ V%j)^ = Trivially, lim„_i.oo ^"^S = i. Thus for any B < \f2 and 

all sufficiently large n, ^-^^ > (v^ ~ V%i)^- It follows that 7^ > (\/ra — \f\if' for all sufficiently 
large n by (69). 

Now, suppose for contradiction that B* < \f2. By Lemma 22 and the fact that B* G (0, 2) by 

i)s2(-i3) 

Lemma 23, we have that d ^ (-g y is a differentiable function of B at B* , and 



Dn2(-B) 



dBDB^_^{-B) 



B=B* 



-1. 



*\2 



)^ since by def. — ; — — — = 

2 ' (-B*) 2 



It follows from definitions and the fact that = \ that 



< 0. 



B=B* 



Thus there exists e > s.t. C{B) < on {B*,B* + e) since C{B*) = 0. It follows that there 
exists B' G (0, \/2) such that C(^') < 0. Thus if we define all relevant functions in terms of B' , 
^n((-v/^ — V%i)^) < for all sufficiently large n by Corollary 2. Also, (-y/n — \/\i)^ < 1 for 
sufficiently large n by Lemma 4. Thus by Lemma 3.i, for all sufficiently large n, 7„ < (-v/n— \/A^)^. 
But this is a contradiction, since we have already shown that B' < \pl implies that 7„ > (-v/ra— \/A^)^ 
for all sufficiently large n. □ 
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Proof of Theorem 12. In light of Lemma 23 and Lemma 24, all that remains to be shown 
is that C{B) < on {B*,2). In light of Lemma 24, it would suffice to prove that C,{B) is strictly 
decreasing on (^/2, 2). But for B G (^2, 2), 



d d Db^{-B) ^ ^ 

——C(B) = — - I — — 

dB^^ ' dBDs^ A-B) dB 2 

< — 7j — 1 H — by Lemma 22 
ij^ 2 

= 0. 

This demonstrates that C{B) is strictly decreasing on (\/2, 2), concluding the proof of Theorem 
12. □ 

Proof of Proposition 1. That B* G [\/2, 2) follows immediately from Theorem 12. That 
Z^i^oo) < min(l, -^) for B > B* follows immediately from Theorem 12 and Lemma 17. i. □ 

7. Proof of Theorem 5, Corollary 1, and Theorem 6. 

Proof of Theorem 5. First, suppose < B < B*. Then i? < 2 by Theorem 12, and thus 
{^/rl — < 1 for all sufficiently large n by Lemma 4. Ci^) > by Theorem 12, and thus since 

B < 2, 4'n((\/^~ V^^)^) > for all sufficiently large n by Corollary 2. It follows from Lemma 3.ii 
that 7„ = {\/ri — for all sufficiently large n. That lim„_i.oo 7n = then follows from Lemma 

4. 

Now, suppose B = B*. Since i?* < 2 by Theorem 12, it follows from Lemma 4 that {^/n — 
y/X^)'^ < 1 for all sufficiently large n. By Lemma 3(i and ii), for all sufficiently large n, either 
^n{{\^ — V%i)^) < 0) which case 7„ = is the unique zero of in the interval 

or 7„ = {^/n - a/A^)^- Let {rii} denote the subsequence of {n} for which 
^n((\/^ ~ < 0. If {ui} is a finite set, then trivially 7„ = (-^/n — \/A^)^ for all sufficiently 

large n. That lim„_>oo 7n = 4 then follows from Lemma 4. 

Alternatively, suppose {rii} is an infinite set. Suppose for contradiction that lim inf j_^oo ^ 

Note that by Lemma 17.n, Z+{^oo) = and ^ooix) > on [0, ■^^^). It follows that there 

exists e > such that < liminfj^^oo ^^(^nj+e < -^^^ < 1 and ^'oo(liminf j^oo ^^(^nj+e) > 0. 
It follows from Lemma 4 and Theorem 10 that for all sufficiently large i, < liminfj^oo ■^"''(^nj + 
e < (\/^~ y/ -^nJ^ < 1 and . (lim inf j_>oo -^^^C^nJ + e) > 0. By the monotonicity of (from 
Lemma l.iv), it follows that for all sufficiently large i, ^'„.(x) > on (— 00, liminfj_>.oo 2'"'"(^'„.) + e). 
But by the definition of liminf, there exists an infinite strictly increasing sequence of integers {n'-} 
s.t. Z+{^ni_) < liminfj^oo ^"^(^nj + e for ah i. Thus for all sufficiently large i, (Z+(^'„/)) > 0. 
This is a contradiction, since ^n'XZ~^i^n')) — for all i. This proves that liminfj_!>oo ■^^'''(^ni) ^ 

4 ■ 

Note that (y^ — A„-)^ < 1 for all sufficiently large i by Lemma 4 and Proposition 1. It 
follows from Lemma 3.i that for all sufficiently large i, Z+(^„.) < (y^ — y^A„.)^, and thus 

limsupj_j.oc Z+(^'„.) < ^ J by Lemma 4. Combining the above, we find that limj_^oo ^ (^nj = 

^ ^' . Thus 7n alternates between two sequences (that of the (yn— VAi) and that of the Z"*"(^„.)), 

both of which converge to ^ . ' -, proving that lim^__^Q(2) — — 4 — — Z (^00 )• 

Now, suppose B G {B*,2). Then (y^ — y/X^)'^ < 1 for all sufficiently large n by Lemma 4. 
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CiB) < by Theorem 12, and thus since B < 2, '^n{{^/n — \/A^)^) < for ah sufficiently large 
n by Corollary 2. It follows from Lemma 3.i that 7„ = Z^{^n) for all sufficiently large n. That 
lim^^oo In = Z~^i^oo) then follows from Theorem 11. 

Now, suppose B > 2. Then {^/n — \/A^)^ > 1 for all sufficiently large n by Lemma 4. It follows 
from Lemma 3.iii that 7„ = Z~^(^n) for all sufficiently large n. That lim„_s>oo 7n = ^^(^oo) then 
follows from Theorem 11. This treats all cases, and we conclude that lim„_^oo 7™ = 7b exists for all 
S > 0. We also conclude that for < 5 < B*, 75 = ^; for B > B*, jb = ^+(^oo)- □ 

Proof of Corollary 1. Suppose for contradiction that liminf„_5.oo v^(l - Pn) < B*. Then 
there exists e > and an infinite strictly increasing sequence of integers {rii} s.t. ^/ni{l — p*.) < 
B* — e, and thus p* . > 1 — "^7^- Consider the sequence {Zi} of continuous time markov chains, in 
which Zi is an M/M/rii queueing system with Xn- = nj — (B* — e)y/ni, (i = 1. Since B* — e < B* , ii 
follows from Theorem 12 that C,{B* — e) > and B* — e < 2. Now, let us define all relevant functions 
(i.e. ^mix)) w.r.t. B* — e. It then follows from Corollary 2 and Lemma 4 that for all sufficiently 
large i, {{^i- \/\^)^) > 0, and {^/n'i- \A^)^ < 1- Thus by Lemma 3.ii, 7^- = {^/r^i- y^A^)^ 
for all sufficiently large i. However, note that by assumption = 1 — < foi^ ^ ^- But this 



is a contradiction, since by Theorem 4, ^ < p* . implies that the spectral gap 7„- of Zi is strictly 

less than — y^A„.)^. 

Alternatively, suppose for contradiction that \^^ms\^^p^^^^/n{l — p*) > B* . Then there exists 
e G (0,2 — S*) and an infinite strictly increasing sequence of integers {ui} s.t. y^(l — p*.) > B* + e, 
and thus p* . < 1 — ^^j^ ■ Consider the sequence {Zi} of continuous time Markov chains, in which Zi 
is an M/M/rii queueing system with A„j = Ui — {B* + €)y/ni, jj, = 1. Since B* + e £ {B*,2), it follows 
from Theorem 12 that C{B* + e) < 0. Now, let us define all relevant functions w.r.t. B* + e. It then 
follows from Corollary 2 and Lemma 4 that for all sufficiently large i, {{y/nl— y^A„.)^) < 0, and 
(\/^ ~ \/ ^rii)"^ < 1- Thus by Lemma 3.i, 7„. = Z+ (*&„.) < (y^ — A„J^ for all sufficiently large 
2. However, note that by assumption = 1 — "^/i^ > pjj. for all i. But this is a contradiction, 



since by Theorem 4, ^ > p* . implies that the spectral gap 7„. of is equal to (y^ — y^ A„.)^. 
It follows that liminf„_>oo y^(l — Pn) = l™sup„_^oo ~ Pn) = B* . □ 

We now prove a series of lemmas, which will collectively prove Theorem 6. To make the depen- 
dence on B clear, let 00(^,2;) denote ^00(2:) defined w.r.t. B. 

Lemma 25. For each B > there exists > s.t. < B' — B < eb implies 75/ > 7^- 

Proof. Suppose < B < B* . Then 7b = ^ by Theorem 5, and thus trivially 75 is strictly 
increasing on (0,B*). 

Suppose B > B* . By Proposition 1 and Theorem 5, 

(70) For all S > B*, < 7B < min(l, — ). Also, < 75. = < 1. 

It follows from (70), Lemma 18. ii. Theorem 5, and the definitions of ^'oo(a;) and jb that for all 
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B > B*, 



y=B 



-IB - (- 



d^b{-b) 



-B) 



2 



A^B + iB'- 47b)^)) -b(^-^{B + {B^- 47b)^)) 
+[1{B + {B^ - 47b)^)) {]^{B - (B^ - 47s)^) j 



-7B + 



(71) 



0. 



It follows from (70) that for all B > B* and e > 0, 
(72) 



'B + e) + {{B + ey-AjB)^ 



5 + (5^-47^)5) > |. 



Thus in light of the definition of ^c>o{x), it follows from (71), (72), and the definition of derivative 
that for B > B* , there exists es > such that B' £ {B,B + es) implies ^oo(-B',7b) > 0. The 
lemma then follows from Lemma 17.i and Theorem 5. □ 



Lemma 26. 75 is continuous on (0, 00). 



Proof. First, suppose B G (0, i?*]. Then 7_b = by Theorem 5, from which it follows trivially 
that 7b is continuous on (0,5*) and left-continuous at B* . 

Now, suppose B > B*. It follows from (70) that there exists 5b G (0,min(l — 75, — ^b,^b)) 

and 5o,B e {0,mm{B,B - B*)) s.t. jb + Sb < min (l, i£z^l^) and B - Sq^b > B* . It follows 
from Lemma 16 that \B' — B\ < 5q^b implies *I'oo(-B', j;) is a differentiable function of x on (7^ — 

5b,1b + 5b) , and for x G hB - 5b,'Jb + • 



dx 



^oo{B\x) < 



< 
< 
< 



(73) 



i{Br-4x)-' 

iB + do,B) 

(B + By'^ since (5o,b < B 
-{4:B^)-\ 



since x > 

-2 



Dx(-y) 



By (70) and the continuity of — i{-y) 
there exists > s.t. \B' — B\ < 5i^B,e implies 



(as a function of y) ensured by Lemma 18.ii, for all e > 



(74) 



D^si-B) 



D 



IB ^ 



-B') 



D. 



7B-1 



i-B) 



D. 



-B' 



< 



Note that for x > 0, ^{B + {B"^ — 4x)2) is a continuous function of B on (2x2, 00), and a strictly 
increasing function of B on [2x2, 00). It follows from (70) that for all e > there exists 52,B,e > 
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s.t. \B' - B\ < 52,B,e implies B' > 2j?,, and 



(75) 



B' 

e 



^{B + {B^ - 47b)^) -\{B' + {{B'f - 47b)^ 



< 



8^2- 



From the definition of ^'oo(a^), (74), (75), the triangle inequality, and the fact that ^'oo(-B,7b) = 0, 
for all e > there exists (^3,5, e > s.t. \B' — B\ < 5^^B,e implies 

(76) \^oo{B',^b)\< ' 



4^2- 

It follows from (73), (76), and the construction oids that for all e G (0, (5b), \B' —B\ < min(5o,B) <^3,B,e^ 
implies ^oo(-B',x) has a zero in (7^ — e,7_B + e). Thus since (7s — e, 7_b + e) is itself contained in 
the interval (0,min(l, ■^)) by construction of 6b, it follows from Lemma 17. i and Theorem 12 that 
this zero must be Z"^(^'oo) (defined w.r.t. B'). Thus by Theorem 5, \B' — B\ < min((5o,_B, (53^B,e) 
implies |7s — jb'\ < £• This concludes the proof that 7^ is continuous on {B*,oo). 

We now prove that 75 is right-continuous at B*. By Lemma 25, there exists 6 > s.t. B' S 
{B*,B*+6) implies 7^/ > 75.. It follows from (70) that for all B' G {B*,B*+6), < 7ij < 

the right-continuity of 7^ at B* then follows from the continuity of Combining the above, we 
find that jb is continuous on (0,oo). □ 

Corollary 3. 7^ is strictly increasing on (0, 00). 

Proof. Suppose < i?i < B2- Let 7max = ^'^Pb£[Bi,B2]^b- Since 7b is continuous and non- 
negative on [Bi, B2] we know that < 7b < 00, and 7max is attained by at least one B G [Bi,B2]. 
Now, suppose for contradiction that there exists some .Bmax S [^1,^2) attaining this maximum. 
Then by Lemma 25, there exists cb^^^ > s.t. B' e (-Bmax, -Bmax + eSmax) implies B' G [-61,-62] 
and 7b' > TSmax = 7max- But this is a contradiction, since by construction 7max = supBg[B^_B2] Jb- 
It follows that the unique B £ [Bi, B2] at which 7^ is maximized is B2, and thus 7^3 > 7Bi- D 



Lemma 27. B > 10 implies < 1 - e-^^^,1 - e-^-8iog(3B+7) ^ min(l, ^), 'I'oo(l 
^^) > 0, and ^00(1 - e"^-8i°s(3B+7)) < q. 



Proof. We proceed by bounding the terms that appear in the representation for ^ooix) when 
one expands the parabolic cylinder functions using their integral representations, as in Lemma 15. 
Suppose B > 10, and .5 < x < 1. Then 



e ' 2 ' ^dz 







< / z ^dz + / e 2 2;(i2; since e 2 < l ^ and z < z for 2; > -B + 1 

< / '-^dz + 6 2 zd^; since 0<z-B<{z- By for z > -B + 1 
Jo Jb+1 

(5 + 1)2-^ 

- ^ ^ +2e-2(3 + B) 



2 - X 

r^+i , (5 + 1)2-^ /-"^ (z-B) _i 

since / z ^dz = - , and / e 2 zdz = 2e 2(3 + 1?) 

7o 2 - X 7b+i 

(77) < (5 + 1)2 + 6 + 25 < {3B + 7f since < x < 1. 
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Similarly, 

e 2 z ^dz > e ^ dz hy non-negativity, and since z ^ > 1 for z > B 

Jb 

R 

(78) = J — > 1 by standard integration techniques. 

Likewise, 

e 22; ^(i^; = e 2 z ^dz + / e 22: ^rf^; 

(B-l)^ Z"-*- _ /""^ (B-z)2 _ (_B-z)2 

< e 2 / a: ^(^2; + / e 22 ^(^2; since i? > 10 implies e 2 inc. on [0, 1] 

< e 2 / 2 ^dz + e 2 dz since z ^<lfor 0<x<l,z>l 

<e 2 z ^dz + e 2 dz by non- negativity 

^0 J~oo 

(B-l)^ (B-l)2 



e 2 . — e 2 

(79) = |-v27r < 1-3 by standard integration techniques. 

1 — X I — X 



Also, 



^ _ ,"1 (B-z)2 

e 22 ^dz > e 2 z ^dz 



-Sl _ . (B-z)2 

> e 2 z ^dz since S > 10 implies e 2 increasing on [0, 1] 

JO 



e 2 

(80) = f . 

1 — X 

It follows from (77) - (80) and Lemma 15 that 

_ (13-1)2 

/ N /C 2 1 Dt(—B) , sO b2 ^ 

(81) (-rTir+^) Sof7(=i) + «^P^ + ^)'=^(i-^)- 

Note that .5 < x < 1 and > 10 implies x < and 

(82) ^((^2 -4x)3 -5) = -2x{{B^ -Ax)^ + B)~^ . 

It follows that the left-hand side of (82) is at least —2B~^, since x < 1 and (i?^ — 4x)^ + B > B. 
Also, the left-hand side of (82) is at most — .5i?~^, since x > .5 and {B"^ — 4x)2 + B < 2B for 
> 10 and x < 1. Combining with (81) and the definition of ^'oo(a^), we find that 



(83) {- ^+3) ^ -2^-^ < ^oo(x) < (35 + 7)^6^(1 -x) - .55-^ 

1 — X 
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Since B>W implies 1 - e-^-8i°s(3S+7) g (.5, l), it follows from (83) that 

_ -81og(3B+7)) < ^ ^^2gf g-f -81og(3B+7) _ ^55-1 

= (35 + 7)"^ - .5B-^ 

(84) < since {3B + 7)"^ < B'^ , B'^ - .bB'^ < for B > 10. 
Similarly, 

(85) ^-00(1 - e"^^^) > ^ - 25"^ > for S > 10. 

The lemma follows by combining (84) and (85). □ 

Proof of Theorem 6. That 7^ is a continuous, strictly increasing function of B follows im- 
mediately from Lemma 26 and Corollary 3. That lim^-^ooTs = Ij and lims^oo ^"^^^^^'^^^ = — i 
both follow from Theorem 5, Lemma 27, and Lemma 17. i. □ 

8. Explicit Bounds on the Distance to Stationarity. In this section we prove Theorem 7 
by demonstrating some novel bounds on the {Qn,kix)} polynomials, and combining with Theorem 
9 and Corollary 1. We begin by bounding \Qn,n{x)\, \Qn,n-i{x)\, and \Qn,n{x) - Qn,n-i{x)\ in terms 
of Cn{x) and hn{x). Let hn{x) = 2nbn{x)~^ . 

Lemma 28. x e [{^/n - \/A^)^, {y/n + \/A^)^) implies 

(i) \Qn,n{x)\Cn{xY'^i < Kix), 
(a) \Qn,n-l{x)\Cn{x)~2 < /l„(x), 

(Hi) \Qn,n{^) - Qn,n-l(a:)|c„(x)"2 < (|)2/i^(x). 

Proof. Expanding Qn,n+i{x) using the recursive definition 12 for the Qn,k{x) polynomials, 

Cji(^) — Qn,n{''^) Qn,n—l{x^ - Qn,n{x) -r (5n,n— 1(^)) 

(86) = Oi^nix) - ^" -Qn,n{x)Qn,n-l{x) + ^Qi^ri-l{x)- 

We now prove assertion i. First, note that 6„(x) > and c„(x) = Qn^ni^) ~ Qn,n-i{x)Qn,n+i{x) > 
by (11) and (13), ensuring that Cn{x)'^ and 6„(x)~^ are well-defined. If Qn,n{x) = 0, then 

\Qn,nix)\Cnix)~2 = < /l„(x). If Qn,nix) / then 

q2 (^j.'jf, {x)~^ (l ■^n ~^ ''T' ~ X Qn,n—l{x) ^ Tl ^Qn,n—l{x) ^2"^-^ 

AjT, Qn,n{x^ Xji Qn,n(,x) 

(87) < sup(l-^^- z + —z^) . 

z An An 

By elementary calculus, 01(2:) = 1 — Aa+Hii^^; -|- j^z'^ is a convex function of z minimized at 
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z = ^ , and 



2n ~ A^^ 2n A^^" 2n 

_ ^ (An + n-x)^ ^ (An + n-x)^ 

2A„n 4A„n 
_ 4Anra - (An + n - x)^ _ fe^(x) 
4A„n 4A„n 

Combining with (87), taking the square root of both sides, and observing that 2\/Xnn < 2n com- 
pletes the proof of i. 

We now prove assertion ii. liQn^n-iix) = 0, then |(5n,n-i(a^)|cn(a^) ^ = < hnix).IfQn,n-iix) / 
then 

q2 (x)c (x)~^ (( '^"'"^'^•^ ^2 Xn -\- Tl — X Qn,n{x) ^ Tl 

' Qn,n—l{x) An Qn,n—l{x^ Xn 

fQQ\ ^ (2 Xn + n-x n .-1 

(88) < sup (z^ z + —) . 

z Xn Xn 

Letting 02 (-z) = z"^ — + y-, we find by elementary calculus that 02 (-2) is a convex function 

of z minimized at z = -'^"+"~^ and 

^A„ + n — x^ ^An + n — x^2 A„ + n — xA^ + n — x^ n 

2Xn 2Xn Xn 2Xn Xn 

_ {Xn + n — x)^ (A„ + n — x)^ ^ n 

4A^ 2Xn Xn 

_ AXnU - {Xn + n- xf _ bl{x) 
4A2 - 

Combining with (88), taking the square root of both sides, and observing that 2A„ < 2n completes 
the proof of ii. 

We now prove assertion iii. It is shown in [38] that Qn,n{x) and Qn.n-iix) do not have any 
common zeros. Thus first suppose Qn,n{x) = 0. Then from (86), c„(x) = x"'3nn-i(^)- Thus 

(89) {Qn,n{^) - Qn,n-l{x)Y Cnix)''^ = ^"'"-^^ ^ = 



Furthermore, 



4A„x 4A„x 



b'i{x) 4Ann - (A„ + n - x)2 

4A„x - 4A„n + {Xn + n- x)^ 



4A„(x - re) + A^ - 2A„(x - re) + (x - n)^ 
hl{x) 



- ' + bl{x) ^ 
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Combining (89) and (90) with the fact that ^ < 1, we find that 

(91) (Qn,n{^) - Qn,n-lix))'^Cn{xy^ < ■ 

bi[x) 

Now suppose Qn,n-i{x) = 0. Then from (86), c„(x) = Ql^^{x). Thus 

(92) {Qn,n{^)-Qn,n-l{x)fcn{x)-^ = ^f^^ = ^ ^ ^by(90). 



Now suppose Qn,n{x) / and Qri,n-i{x) / 0. Then 



i^Qn,n{x) Qn,n— 1(^)) C-nix) 



Qn,,n(x) _ -^-^2 

/ Qn,n(x) \2 _ Xn+n — X Qn,n{x) , n_ 

^Qn,n-l{x) ^ A,i Qn.n-l{x) A„ 

(z-l)2 



(93) < ^^P3 A„,+n-.^| n- 

Let 03(2;) = ^ — A^+n-l inf2GM('^^ ~ An+n-x- ^ _^ ^) > from our earlier analysis of a2(^)- 

Thus by elementary calculus 03(2) is a continuously differentiable rational function of z on M, and 

— a3(z) = ^ ^ — 

^ An, An ' 

{z - l)(2z2 - 2 ^"+^-^ 2 + - 2^(2 - 1) + - 1)) 

^ An An 

_ (z - l)((n - A„ + x) + (A„ - n + x)z) 
A„(z2-^4ll^z + ^)2 ■ 

It follows that the zeros of ^03(2;) occur at z = 1 and 2; = Thus by elementary calculus, 

sup^a3(z) must be one of 03 ( 1) , as ( ) , lim^^^oo 03 (^) , lim^-^oo Q3 (^) ■ Trivially, 03(1) = 0, 
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and lim2_j>_oo 03(2) = lim2_!.oo a3{z) = 1. Furthermore, 



Xn-n-x 2x 
aslT : — ) = 03(1 



Xn — n + X Xn — n + X 

\n—n+x 



( 2x \2 
V —n+x ) 



2x \2 



( ^ )^ 

^ X„-n+x ' 



Ax _|_ / 2x \2 Xn+n—x < Xn+n—x 2x i _n_ 

Xn—n+X ' ^Xn~n+X ' Xn Xn—n+xXn Xn 



( 2x \2 
yXn-n+x' 



I _ 4a: I / 2x \2 _ Xn+n-x i Q i 2{n-x) n 2x , _n_ 
Xn-n+x \Xn-n+x) Xn V I X_^-n+X ) Xn Xn 



Ax 

[Xn-n+x)'^ 



-A I Ax _| 3_ _|_ An- Ax 

Xn-n+x {Xn-n+x)'^ A„ ' A„(A„-n+i) 



AxXr 



-4A„(A„ -n + x) + AxXn + 4(n - x)(An - n + x) + 3(An - n + x)^ 

4xAn 

(A„ - n + x) ( - 4A„ + 4(n - x) + 3(A„ - n + x)) + 4xAn 

4xAn 

(A„ - n + x)(-A„ + n - x) + 4xA„ 
4xAr, 4xAr, 



(94) 



Combining the above with (90), we find that s\XY)^az{z) < ^^"^ 



4xA„ - (A„ - n + x)2 b'lix) 



sup 03(2:) < max(a3(l), a3(— ), hm 03(2), hm 03(2;)) 

z Xn — n + X z^oo 

= max(0, -^^^,1,1) 

Combining with (91) - (93) and taking square roots completes the proof of the lemma. □ 

We now prove some simple asymptotic bounds that are necessary to proceed with our analysis. 
Recall that g„ = 1 + n^^. Let r„(x) = qnxn~^, and s„(a) = qnan~2. Then 

Lemma 29. For each a > B > 0, there exists Nq. s.t. n> Na, /c G [n — a^/n — 3, n + a^/n + 2>\, 
and X > implies 



I ^r- I I — I 

I min(A;,n) I ' I A„ I 
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Proof. For any constants a,b > and all sufficiently large n, 
n _ an2 + b 



1 + 



n — an2 — b n — an^ — b 



1 (a + b)n + abn2 



^-,2 

— 1 + an~2 + - 3 

— an2 — bn 

= l + an"3 +0(n"^) 

_ 1 _ 3 

< 1 + an 2 + an 4 for all sufficiently large n 

(96) = l + Sn{a). 

Similarly, for any constants a,b > and all sufficiently large n, 

(n-oni-6)"^ < (1 + s„(a))n~Hy (96) 
= n + g„.a« 2 

< n^"*^ + n~3 for all sufficiently large n, since qnan~^ = 0{n~~^) 

(97) = qnn-\ 

Thus from the assumptions on k, the definition of A^, and the assumption that a > i?, we find that 
for all sufficiently large n, 

(98) min(A;,n)~"^ < qnn~^ by taking a = a,6 = 3 in (97); 

(99) < qnn~^ by taking a = a, 6 = 3 in (97) since A^^ < (n — ana — 3)~^; 

(100) r- < 1 + Sn(a) since A^ < n, and applying (96) with a = a, 6 = 4; 

mm[k — 1, n) 

n 

(101) T— < 1 + Snia) by applying (96) with a = a,b = 0. 

An 

Assertions i. and ii. follow immediately from (98)-(101) and the assumed positivity of x. □ 

We now bound \Qn,k+iix)\ and \Qn,k+i{x)-Qn,kix)\ in terms of\Qn,k{x)\ and \Qn,kix)-Qn,k-iix)\ 
and visa-versa, which will later enable us to bound \Qn,kix)\ for k = ni: 0{y/n) using Lemma 28. 
More precisely, we prove 

Lemma 30. For each a > B > 0, there exists s.t. n> Na, k G [n — a^Jn — 3, n + a-^/n + 3], 
and x > implies 

(i) \Qn,k+l{x)\ < {l + rn{x)){l + Sn{a))[\Qn,k{x)\ + \Qn,k{x) - Qn,k~l{x)\) , 
(a) \Qn,k+l{x) - Qn,k{x)\ < {I + rn{x)){l + Sn{a)){rn{x)\Qa,k{x)\ + \Qn,k{x) - Qn,k-l{x)\) , 
(Hi) \Qn,k-l{x)\ < {l + rnix)){l + Snia)){\Qn,k{x)\ + \Qn,kix) - Qn,k+l{x)\), 

(iv) \Qn,k-l{x) - Qn,k{x)\ < (l + rn{x)) (l + ^^(q:)) {rn{x)\Qn,k{x)\ + \Qn,k{x) - Qn,k+l{x)\) ■ 



37 



Proof. Let us choose Na to satisfy the conditions of Lemma 29 (the existence of such an 
foUows from the same lemma), and suppose n > N^- We first prove assertion i. Recah that by for 
k >1, \Qn,k+i{x)\ is equivalent to 

X mm(k,n) , , min(/c,n)^ , 

1(1 - 1- + ^^)Qn,kix) ^^Qn,fc-l(x)| 

I/-, ^ / N min(A;,n).„ , , „ , 

= 1(1- —)Qn,k{x) + (Q„,fc(x) - Qn,k^l{x)) I 

< |1 - ■T-\\Qn,k{x)\ + I — ^ — \\Q7i,k{x) - Qn,k-i{x)\ by the triangle inequality 

< (l + rn{x))\Qn,k{x)\ + (l + Sn{a))\Qn,k{x) - Qn,k~i{x)\ by Lemma 29. i and Lemma 29. ii, 
from which assertion i. follows, since by the strict positivity of r„(x) and s„(a) we have 

(102) (1 + r„(x)) < (1 + r„(x)) (l + s„(a)) , and (l + s„(a)) < (l + r„(x)) (l + s„(a)) . 

We now prove assertion ii. Rearranging definition (12), we find that for A; > 1, \Qn,k+i{x) — 
Qn,k{x)\ is equivalent to 

. X ^ , . min(A;,n) , . 

I - -^Qn,k\x) + 7 [Qn,k{x) - Qn,k-l{.x))\ 

< I - ^ll<3n,fc(a::)| + I — WQnM^) - Qn,k-i{x)\ by the triangle inequality 

< rn{x)\Qn,kix)\ + (l + Snia)) \Qn,k{x) - Qn,k-iix)\ by Lemma 29.i and Lemma 29.ii, 
from which assertion ii. follows by (102) and the similar fact that 

(103) rn{x) < (1 + r„(x)) (l + s„(a))r„(2;). 

We now prove assertion iii. Rearranging definition (12), we find that for k > I, \Qn^k-iix)\ is 
equivalent to 

I ■ x ((l - ^ + T )Qn,kW -Qn,k+l[x)}\ 

mill( /C, Tl I 

= K^^^T V ^1 T + l)Q™,fc(2;) r^T 

min(«;,n) mm(A;,nj mm(A;,nj 

X \ 
= 1(1 ^)Qn,k{x) H — r{Qn,k{^) - Qn,k+l{x))\ 

mm[k,n) min(A;,nj 

< 1 1 — \\Qn,k{x)\ + \ . " — r \\Qn,k{^) - Qn,k+i (x) \ by the triangle inequality 

mm(A;,n) mm(A;,nj 

< (l + rn{x))\Qn,kix)\ + (l + Snia)) \Qn,kix) - Qn,k+iix)\ by Lemma 29.i and Lemma 29.ii, 

from which assertion iii. follows by (102). 

We now prove assertion iv. Rearranging definition (12), we find that for A; > 1, \Qn,k~iix) — 
Qn,k{x)\ is equivalent to 

A„ , X mm(k,n) 

((1 - T- H ^ )Qn,k{^) - Qn,k+l{x)) - Qn,k{^)\ 



min(/c, n) A„ A 
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-)Qn,k{x) r^T ^Qn,k+l{x) 



'min(A:,n) min(/c,n) ' min(A;,n) 

-Qn,k{^) \{Qn,k{^) - Qn,k+l{x))\ 



min(A;,n) ' min(A;,n) 

< I — — - \\Qn,k{x)\ + \ . " — r \\Qn,k{x) - Qn,k+i {x) \ by the triangle inequality 

mm[K,n) min(K,nj 

< rn{x)\Qn,k{x)\ + (l + Sn{a))\Qn,k{x) — Qn,k+i{x)\ by Lemma 29. i and Lemma 29. ii, 
from which assertion iv. follows by (102) and (103). □ 

We now use Lemma 30 to inductively bound \Qn,k{x)\ for values of k that are 0(n2) away from n. 

Lemma 31. For each a > B > 0, there exists Na s.t. n > Na, and x G ((a/^ — iV^ + 

-y/A^)^) implies that 

(i) For /c € [n, n + a^/n + 3], 

\ k—n 

(l + r„(x)) (l + Sn{a)) (l + r„(x)2) j c„(x) 2 /i„(a;), 

(h) \Qn,k{x) - Qn,k~i{x)\ < ((1 + rn{x)) (l + s„(a)) (1 + r„(x)^)) r„(x)5c„(x)3/i„(x). 
('iij For € [n — a^/n — 3, n — 1], 

(1 \ 1 
(1 + r„(x)) (1 + Sn{a)) (1 + r„(x)2) j c„(x) 2 /i„(a;), 

/ 1 \ 1 1 

(b) \Qn,k{x) - Qn,k+i{x)\ <[{}+ r„(x)) (1 + s„(q)) (1 + r„(x)2) j r„(x)2c„(2;)2/i„(x). 



Proof. Let us choose Na to satisfy the conditions of Lemma 30 (the existence of such an 
follows from the same lemma), and suppose n > N^- We first prove assertion i. We proceed by 
induction on ia. and ib. simultaneously. The base case k = n follows immediately from Lemma 
28. i and Lemma 28.iii. Now, suppose the induction is true for some k £ [n,n + an^ + 2]. Then by 
Lemma 30. i, 

\Qn,k+lix)\ < (1 + r„(x)) (1 + S„(a)) {\Qn,k{^)\ + \Qn,k{^) " Qn,k-l{x)\) 

< (l + r„(x))(l + s„(a)) ^((l+r„(x))(l + s„(a))(l + r„(x)5))^ \(x)5/i„(2;) 

+ (^(1 + r„(x)) (1 + s„(a)) (1 + r„(x)5)^ r„(a;)3c„(x)5/i„(a;) J by the ind. hyp. 



= ((l + r„(x))(l + s„(a)))''^' "(l + r„(x)5)^ "c„(x)U„(x)(l + r„(x)^) 
(104) = (^(l + r„(x))(l + s„(a))(l + r„(x)2)j c„(x) 2 /i„(x). 
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Similarly, by Lemma 30. ii 

\Qn,k+lix) - Qn,k{x)\ < (1 + (1 + Sn{a)) (r„ (x) | Q^^fc (x) | + |Q„,fc(x) - Qn,k-i{x)\) 

< (1 + r„(x)) (1 + Sn{a)) ^((1 + r„(x)) (l + s„(a)) (l + r„(x)^))^ "r„(x)c„(x)U„(x) 

+ (^(l + r„(x)) (l + s„(a)) (l + r„(x)5)^ r„(x)5c„(x)5/i„(x)^ by the ind. hyp. 

= (^(l + r„(x))(l + Sn(a))) (l + r„(x)^)''""c„(x)5/i„(j;)(r„(x) +r„(x)^) 

/ 1 \ fc+i— " 1 1 

(105) = (^(l + r„(x))(l + s„(a))(l+r„(x)2)j r„(x) 2c„(x) 2 



This concludes the induction, proving assertion i. 

We now prove assertion ii., which follows from a very similar argument. We proceed by induction 
on iia. and iib. simultaneously. The base case k = n — 1 follows immediately from Lemma 28. ii and 
Lemma 28.iii. Now, suppose the induction is true for some k £ [n — an^ — 2, n — 1] . Then by Lemma 
30.iii, 

\Qn,k-l{x)\ < {l + rn{x)){l + Sn{a)){\Qn,k{x)\ + \Qn,kix)-Qn,k+l{x)\) 

< ((l + r„(x))(l + s„(a))(l + r„(x)5))" ^^„(x)U„(x) 
by exactly the same argument used to derive (104). 

Similarly, by Lemma 30. iv, 

\Qn,k-lix) - Qn,kix)\ < (1 + Vnix)) (l + S„(a)) (r„(x) | Q^^fc (x) | + \Qn,k{x) - Qn,k+l{x)\) 

< (^(l + r„(x))(l + s„(a))(l + r„(x)5)) r„(x)5c„(x)U„(x) 
by exactly the same argument used to derive (105). 
This concludes the induction, proving assertion ii. and thus the lemma. □ 

Lemma 32. For each a> B > 0, there exists Na s.t. n > Na, x G ((v^-\/A^)^, (^Al+^A^)^), 
and k €z [n — a^/n — 2, n + a^/n + 2] implies that 

(i) \Qn,k{x)\cn{xY'^2. < /i„(x) exp ((a + n"'l)(a + 3x5)) 

(ii) \Qn,k+i{x) - Qn,k{x)\cn{xy2 < /i„(x)r„(x)2 exp ((a + n"5 )(a + 3x2)) 

Proof. By Lemma 31. ia and Lemma 31. iia, we find that there exists Na s.t. n > Na implies 

\Qn,ki^)\cnixy^ < /i„(x) exp (|n - A;|(r„(x) + s„(a) + r„(x)i)) 

< /i„(x) exp {{an^ + 2)(r„(x) + s„(a) + r„(x)^)) 

= hn{x) exp ^(an2 + 2) (g„xn~^ + qnan~2 + q^x^n'z'j^ 

< /i„(x) exp ^(7„(an2 + 2) (xn^^ + an~2 + 2;2n^2)^ since qn < q-n 

(106) = /in(x) exp ^(7„(a + 2n^2) (xn~2 + Q + x2)^ . 
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By identical reasoning, 



(107) \Qn,k+i{x) - Qn,k{^)\'^nix) 2 < /i„(x)r„(x) 2 exp (a + 2n l){xn 2+0 + ^2' 



Furthermore, note that since 



1 



_ 1 1 

(108) = — < < 1, 

2x2 2?i2 2n2 

it win be the case that 

(109) < X < (\/n + \/A^)^ imphes xn~^ < 2x2. 
Also, note that 

qn{a + 2n~^) = (1 + n~3)(a + 2n~i) 

_1 _i _3 

= a + an 4 + 2n 2 + 2n 4 
= a + 0{n~^) 

(110) < a + n^5 for all sufficiently large n. 

The lemma follows by combining (106), (107), (109), and (110). □ 

We now derive a variant of Theorem 9 for -P"<j(t) (as opposed to -Pj" (t)) that does not simply 
sum over all j + 1 states < j, but instead relies on a 'probability flow' interpretation using the 
Chapman-Kolmogorov equation. This will allow us to bound the distance to stationarity uniformly 
in n for the c.d.f., since we can sidestep having to sum over 'too many' error terms. 

Lemma 33. For ^ > p*, 
(111) 

|J^<,(t)-P^,(00)| < {2TT)-^gn{j)—n'^ / e-^'x-^\Qn,i{x)\\Qn,j+l{x)-Qnj{x)\bn{x)Cn{x)-'dx. 

Proof. Note that 

A;=0 

= /^i(t) - XnP;:,{t) + ^ (min(/c + l,n)J^fc+i(t) + A^/^^iW " + min(A:, n)) J^,(t)) 

k=l 

by the Chapman-Kolmogorov equation 
= mm{j + l,n)P^j_^_^{t) — XnP^jit) by a telescoping sum argument. 

It follows that for ah i,t>0, P^j{oo) = Pl<j{t) + (min(i + l,n)P^^j^^{s) - XnP^j{s))ds. Thus 
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(min(j + l,n)f^«+i(s) - A„/^';,.(s))ds| 
jj"" (min(i + - /^^.,i(oo)) - Xn{P^jis) - P,^,{oo)))ds 



since the steady-state terms cancel by detailed balance 

e~'^''Qn,iix)QnJ+lix)bnix)Cnix)~^dx) 



^ (^min(j + l,n)((27r)-^(7n(i + l)^(A„n)-i 



An((27r)-i5n(j)^(Ann)-^ 



(v^+v^)^ 



(v^-VA;r)2 

e~'^''Qn,i{x)Qn,jix)bn{x)cn{xy^dx) ds 



by Theorem 9 



r ((2vr)-i5n(j)^n-i / 



'Qn,i{x){Qn,j+l{x) - Qn,j{x))bn{x)Cn{x) ^dxjds 



7x1 1 Tl\ 1 Tt! 

since min(j + 1, n)s-„(j + 1) — (A„n)- = Xngnij)Tz^{Kn)~ = gn{j)j^n- 

■^n 



< 



1 r{^/n+^/X7^)^ 



r((2.)-V(4n-| 

(27r)-i5n(i)^n-i / 

by Tonelli's Theorem 

,1 _ r{V^+V>^f 



'\Qn,i{^)\\Qn,j+l{x) - Qn,j{x)\hn{x)Cn{x) ^dx]ds 



e '^''\Qn,i{x)\\Qn,j+l{x) - Qnj{x)\bn{x)Cn{x) ^ ds dx 



(27r) gn{j)—:n 



^|(5n,i(2;)||Qn,i+i(2;) - Q„,j(x)|6„(x)c„(x) ^dx, 



and the lemma follows. □ 

We now prove bounds on a special type of integral that arises in the analysis of Pj^j{t) — P"(oo). 

Lemma 34. For all B,C > there exists Nb,c s.t. n > Nb,c implies that for all c € (0, C] and 
t > 1, 

/ exp ( - xt + cx2)6„(x) ^dx < )2 exp (-Sc + c^(2t) - —t). 

J{V^-Vx^r tXn 2 4 

Proof. By Lemma 4, {y/n — vXi)^ < 2a/A^ < (v^ + VXi)^ for all sufficiently large n. Thus 



(v^-V^)2 



exp ( — xt + cx 2 ) (4A„n — (An + n — x)^) ^ dx 



2\/ A„n 



(112) = / exp(-xt + cxi)(x- (v^n- v^)^) HiVn + ^/^f-x) ^^dx 

(113) +/ exp(-xt + cx2)(x- (\/n- VA^)^) ^ {{Vn + yO^f - x) ^dx by (11). 
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A — — — — 

Let Un = 2y/X~n - - y/K,)'^. Since {{^/n + - x) ^ < {{^/n + ^A^)^ - 2y/X~n) ^ for 

X G ((-y/n — V%i)^i 2\/A^) , the term appearing in (112) is at most 



{{Vn+ VKf -"^VKny^ / exp{- xt + cx2){x- {^/n- ^/Xrlfy^dx 



substituting y = x — {y/n — \/A^)^ 
< (A„ + n) 2 exp - a/A^)^)*) / exp - + cy2 + c(A/n - \/A^))y~^c?y 





111 
n 



since a, 5 > implies (a + 6) 2 < a2 + 6s 

_1 / N r^n ^ ^ 

(114) = (A„ + n) 2exp(^-(V^- VA^)^t + c(^/^- VA^)jy exp (- yt + cy 2)^-2 dy. 
We now show that 

(115) -yt + cy^ < -^yt + c^{2ty^ for all y > 0. 

Indeed, by elementary calculus the concave function a4(y) = {—yt + cy2) — [—^yt+c'^{2t)~^) attains 
its maximum when y = (f )^, and a4((^)^) = 0, from which the desired inequality follows. Plugging 
(115) into (114), we find that the term appearing in (112) is at most 



_1 / X /""n I J 

[Xn + n) 2 exp - ( - VA^)2t + c(V^ - VA^)j exp ( - -yt + c2(2t)-i)y-2dy 



= (A„ + n) exj>l^- {y/^- y%,ft + c{y/^- y%^) + c\2t)-^j exp { - -yt)y-2dy 

_i / N r°° I J 

< {Xn + n) 2exp(^-(V^- VA„)2t + c(V^- VAn) + c^(2t)"^jy exp { - -yt)y-2dy 

= (A„ + n)"^ exp ( - (\/^ - VKft + c{V^ - + c^{2ty^^ (y )5 

since / exp ( yt)y~^dy = ( — )2 by standard integration techniques 

Jo 2 t 



(116)< {^)^exp(-{^-^/X^ft + c{^-^/K) + c^{2t)-^) since A„ + n > 2A„. 

Ant V / 

We now bound the term appearing in (113). Since (x—(y^—-v/\i)^) 2 < (2y/Xn'n—{y/n— -v/A^)^) ^ 
for X G (2\/A„n, (y^ + \/A^)^), the term appearing in (113) is at most 

'2\/A„n 

"(v^+v^) 



(2A/A„n- (Vn- VA„)^) 2 / exp ( - + ca;2 ) ((^ + y^A„) - x) 2(^3; 

J2\/A„n 

1 1 Hv^+v^)^ _i 

< (2vAnn-(A/n-V A„) ) 2 gup exp(-zt+cz2) / {{V^+V^n) -x) ^dx 
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(2vA^ - (Vri - vA^)^) ^ sup exp[-zt + cz^) / y'^dy 

ze\2^/x~^.(^/E+^/x::)^] Jo 



substituting y = {y/n + \/A^)^ — x 
< 2(— )2 sup exp(-zt + cz2 

3A„ - n z^[2y%::^,(V^+y%^y] 



since 2\f\ji — (y/n — \/A^)^ > 3A„ — n, and / y~^dy = 2{Xn + n)2 

Jo 

(117) < 2(-^±^)iexp(c2(2t)-i-(A„n)h) 

by applying (115) to sup exp ( — + czz). 

2:e[2VA„n,(v^+v^)^] 

Furthermore, 

(118) {— )2 < 2 for all sufficiently large n, since lim — = 1. 

3A„ — n n-¥oo 3A„ — n 

Also, 

1 B^c 

c{\/n — \/ Xn) < —Bc-\ J by Lemma 4 

2 4n2 - 2B 

= ^Bc + 0{n~^) 

(119) < —Bc + n~3 for all sufficiently large n. 
In addition, 

(A„n)i > A„ 

1 1 B^ 

(120) > -Xn + An + — for all sufficiently large n, 

I B^ 1 
smce Xri + — = 0{Xr\). 

Combining (116), (117), Lemma 4, and (118) - (120), we find that for all sufficiently large n and 
t > 1, 

/ exp — xt + cx^)bn{x)^^dx 

J{V^~V^)^ 
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< exp (c2(2t)-i) exp [-^t + ^Bc + n's) + 4exp ( - ^A„t - Alt - ^t) 

< exp {c^{2t)-^ ~ ^ ^ '^^^^ (^17^^ + 4exp ( - ^A„t - Ali 



since exp(— i?c + n > 1 



< exp (c2(2t)-i -^t + \bc + n-i) + A{\Kt)-\\lt)-') 

= (— )2exp(c2(2t)-i - —t + -Bc + n-3) 1 + 8(^2 A„t2) 



(121) < (-^)3exp(c^(2t)~^ - ^t + ^Sc + n"3)(l + 8A^^^) since t > 1, tt^ > 1. 

Furthermore, since 1 + 8A^^ = 1 + ©(n^-*^), and exp(n~3 ) = 1 + 0(n~3) by a simple Taylor series 
expansion, it follows that (1 + 8A~^) exp(n~ 3 ) = 1 + 0(n~3). Thus for all sufficiently large n, 
(1 + 8A~^) exp(n~ 3 ) < 1 + n^i = Combining with (121), and noting that in all cases 'for all 
sufficiently large n' can be defined in terms of C only (as opposed to c), completes the proof of the 
lemma. □ 

Proof of Theorem 7. Suppose i? G (0,5*) , 01,02 e M. Let a = max(i?, |ai|, |a2|). Then for 
all sufficiently large n and all t > 1, 

(122) ^\P\n+a,V^,\n+a.^ii) - P\n+a.V-n^^^)\ 

< (27r)~^g„exp (B^)(A„n2)"i / exp ( - xt) |(5n,i(a:^)|cn(x)" 2 |Q„j(x)|c„(x)" 2 6„(x)dx 
by Theorem 9, Corollary 1, and Lemma 5 

< (27r)^^g„exp (i?^)(A„n2)^^ / exp ( — xt)/in(a:;)^ exp (2(a + n"5)(a + 3x2 ))6„(x)(ix 



by applying Lemma 32 to both \Qn,i{x)\cn{x) 2 and \Qn,j{x)\cn{x) 2 

(123) 

= 27r~"^gn exp (i?^) — / exp ( — xt + 2(a + n~5)(a + 3x2 ))6^(x)~^(ix. 

An J{y^-^)2 

We now bound the integral appearing in (123), which is equivalent to 

_i f<^+^)^ _i 1 

exp (2a +2an 5") / exp ( — xt + 6(a + n 5)x2)6n(x) dx 



(124) < g„(^)^ exp (la^ + 2an'^ + W{a + n^i) + 18(q + n~^)^t"^ - 



' V 4 
for all sufficiently large n, by applying Lemma 34 with C = 6(a + 1), c = 6(a + n~5). 
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Also, for all sufficiently large n, 

(125) exp ((2a + 3B)n^5^ = 1 + 0(n~5 ) < 1 + = Qn hy a simple Taylor series expansion, 
and similarly 

(126) exp ((36an"^ + ISn'^t"'^) = 1 + 0(n"i) < 1 + n"3 = g„ since t > 1. 

Combining (123), (125), and (126) we find that for all sufficiently large n, and all t > 1, the integral 
appearing in (123) is at most 

(127) qU^)^ exp (20^ + Wa + ISa^t'^ - ^t). 
Combining (123) and (127) we find that (122) is at most 

27r"^g„ exp {B^) ^ ( g^(^)^ exp (2a^ + 3Ba + ISo^t"^ - 

= 2(7rt)~^g4(iL)| gxp (b^ + 20^ + 3Ba + ISa^t"^ - 

A„ V 4 

(128) < 2(7rt)~^g4( " )f gxp('2(a + 5)2 + 18a2t"^-^t). 
Furthermore, 

1 

13tI 2 1 

(129) ~ -'^ ~ = 1 + 0(n~2) < g„ for all sufficiently large n. 

The first part of Theorem 7 follows from (128) and (129), by noting that 

(130) g^-^ = l + 0(n-4) < 1 + n 8 for all sufficiently large n. 
Similarly, for all sufficiently large n and all t > 1, 

(131) |f^<,.(t)-P^,.(oo)| 
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< (27r) ^g„exp(B^)n W exp ( - xt)a: ^|(5n,i(3;)|c„(x) ^\Qn,j+i{x) - Qn,j{x)\cn{x) ■^hn{x)d2 



by Lemma 33, Corollary 1, and Lemma 5 



r(V^+V>^Y j_ _i 1 _ 

< (27r)~^g„exp / /i„(x)^r„(x) 2 exp ( - xt + 2(a + n~5 )(a + 3x2 

by applying Lemma 32 to both \Qn,i{x)\cn{x)~"2 and \Qn,j+i{x) - Qn,j{x)\cn{x)~"2 
= 27r g^exp(i? )n2 / exp ( - + 2(q; + n 5)(a + 3x2))x 2&„(x) dx 

3 1 /■{v^+V^)2 ^ J 

< 27r-^g^ exp (^2)^2(25"^) / exp ( - xt + 2(a + n's )(a + 3x2 ))6„(x)"^(ix 



'(v^-V^)2 

since x~2 < Uy/^ _ y'^)^) ^ = {^/n — ^/Xn)^'^ < 2B~^ by Lemma 4 



< 4(S7r)-^g|exp(s2)n3^g^(^)lexp(^2a2 + 3Sa + 18a2t-i-:^t)j By (127) 

= 4qf^^{B^TTty^{^)^ exp (^^ + 2a^ + 3Sa + ISa^r^ - ^t) . 

A„ V 4 / 

The second part of Theorem 7 follows by combining the above with (129) and (130). □ 

8.1. Comparison to other bounds from the literature. In this subsection we compare the bounds 
from Theorem 7 to two other explicit bounds given in the literature [46], [7]. In both cases we will 
prove that the bounds from the literature (applied to |-P^<„(t) — P^^{oo)\ for < B < B*) grow 
with n, rendering them impractical in the H-W regime. We begin with the bounds given in [46], 
which translate to the statement that for each B G (0,B*) there exists Nb s.t. n > Nb implies 
that for all t > 0, 
(132) 

CO 

|J'n",<nW - ^<n(00)| < 4(n- 1) ( ^ ((-^)^ - l)i^"(00) + ((-^^ 



i=l 



Note that since lim„_j.oo ^^i^ — 0? the exponential rate of convergence demonstrated by (132) 
goes to zero as n — )■ 00, rendering the bound in [46] ineffective. We now examine the bounds given 
in [7], which translate to the statement that for each B S (0,B*) there exists Nb s.t. n > Nb 
implies that for all i > 0, 

(133) |P„"<„(t) - P^„(oo)| < (P^(oo)-i - 

It is well-known (see [22]) that liminf„_).oo ^^"^'^\ > 0. It follows that the prefactor demou- 

nt 

strated by (133) diverges as n — )• 00, rendering the bound in [7] ineffective. 

It should be noted that although the bounds given in [46] and [7] are ineffective for the particu- 
lar events of interest in this paper, both bounds hold in much greater generality, and thus remain 
interesting and applicable in a variety of other settings. 

9. Conclusion and Open Questions. In this paper we studied the rate of convergence to 
stationarity of the M/M/n queue in the H-W Regime. We explicitly computed the limit of the 
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exponential rate of convergence to stationarity (spectral gap). We proved that there is an inter- 
esting phase transition in the system's behavior, occuring when the excess parameter B reaches 
B* ~ 1.85772. For B < B* , the exponential rate of convergence is above B* it is the solution to 
an equation involving the parabolic cylinder functions. We showed that this transition asymptoti- 
cally characterizes a phenomenon previously observed to occur for fixed n, unifying and simplifying 
several earlier lines of work. We also explicitly bounded the prefactor governing the rate of con- 
vergence uniformly over n (for B < B*), and used our bounds to offer managerial insight into the 
operation of a call center. 

This work leaves several interesting directions for future research. There are many open questions 
related to the interaction between weak convergence and convergence to stationarity. Although our 
results and those of [30] show that for the M/M/n queue in the H-W regime there is an 'interchange 
of limits' in this regard, it is unknown to what extent such an interchange must hold in general. It 
would also be interesting to prove that a phase transition occurs in other models, such as queues 
with abandonments and non-Markovian service distributions. 
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